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1.  INTRODUCTION 


Sharp  gradients  in  the  wind  encountered  by  an  aircraft  on 
takeoff  or  on  final  approach  have  caused  several  serious,  in  some 
cases  fatal,  accidents  to  occur.  Of  the  three  components  of  wind 
as  seen  from  the  aircraft,  namely  headwinds,  crosswinds,  and  down- 
drafts,  gradients  in  headwinds  appear  to  be  the  single  most 
dangerous  condition.  Sudden  crosswinds  have  contributed  to  acci- 
dents by  distracting  the  pilot  from  vertical  position  instruments 
(e.g.,  NTSB,  1973c),  but  have  not  been  identified  as  the  cause  of 
a serious  accident.  Sudden  changes  of  headwind  are  critical  be- 
cause of  the  relation  between  airspeed  and  lift,  and  are  doubly 
insidious  because  of  the  lack  of  distinctive  accelerative  cues 
felt  by  pilot  and  crew.  In  several  accidents  pilots  were  seeking 
to  acquire  visual  contact  with  runway  indicators  at  the  time  that 
instruments  were  indicating  descent  below  the  glideslope  and  ex- 
cessive sink  rates;  even  after  visual  contact  was  obtained, 
neither  visual  nor  accelerative  cues  wer^  distinctive  enough  to 
depict  the  gravity  of  the  situation  in  time  (NTSB,  1974b, c, 

1967a, c) . 

There  are  three  meteorological  conditions  which  are  frequently 
accompanied  by  regions  of  severe  wind  shear  (Langweil,  1976). 

a.  Boundary  layer  shears -- inversions  and  jets;  character- 
ized by  horizontal  strata  existing  for  several  miles, 
but  having  a vertical  profile  with  significant  shears 
(vertical  shear);  while  these  can  be  quite  abrupt  and 
reach  significant  values  (e.g.,  20  kphf  (see 
glossary)  or  0.3  m/sec/m),  they  have  yet  to  be 
positively  identified  with  a serious  accident. 

b.  Frontal  shears --cold  and  warm;  fronts  causing  hazardous 
shears  are  characterized  by  abrupt  wind  shift  zones 
having  considerably  different  wind  directions  on  either 
side.  Warm  fronts  have  shallow  slopes  compared  with 
approach  and  departure  paths,  and  thus  present  a vertical 
shear  region  to  an  ascending  or  descending  aircraft;  cold 
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fronts  have  steep  slopes  in  the  atmospheric  boundary 
layer,  so  that  even  aircraft  flying  at  constant  altitude 
may  encounter  shear  zones  (horizontal  shear)  , The 
Iberia  Airlines  accident  (NTSB,  1974c)  appears  to  have 
been  caused  by  a warm  front  shear. 

c.  Thunderstorms - -gust  fronts  and  downburst  cells;  char- 
acterized by  localized  zones  of  rapidly  varying  wind 
currents,  which  an  aircraft  experiences  as  downdrafts  and 
sudden  changes  in  headwind  and  crosswind.  These  have 
been  identified  as  a prime  contributing  factor  in  several 
accidents  (NTSB,  1974,  1976). 

The  need  exists  for  a sensor  which  can  remotely  detect, 
identify,  and  track  such  wind  shears  around  the  airport  and  pro- 
vide safety  assurance  to  both  departing  and  approaching  aircraft 
up  to  1500  feet  (450m)  or  more  in  altitude.  Point  sensors,  such 
as  anemometers  and  pressure  jump  sensors,  and  vertical  looking 
probes,  such  as  acoustic  sodars,  FM/CW  radars,  and  CW  laser  radars, 
must  be  placed  off  the  airport  surface  if  they  are  to  detect  gust 
fronts  or  downburst  cells.  There  are  two  technologies  which  can  be 
placed  on  the  airport  surface  and  provide  the  potential  of  re- 
motely sensing  shear  events  beyond  the  airport  boundaries:  pulsed 
Doppler  radar  operating  off  returns  either  from  precipitation  or 
index  of  refraction  variations,  and  pulsed  Doppler  laser  xadar, 
operating  off  returns  from  precipitation  or  particulate  aerosols. 
Both  are  promising  techniques.  Lasers  have  now  reached  a state  of 
maturity  that  makes  the  feasibility  of  a combined  vertical  probe 
and  remote  sensor  possible;  it  is  given  primary  emphasis  in  the 
later  sections. 

Section  2 discusses  the  meteorology  of  wind  shear  events. 
Section  3 develops  the  major  requirements  that  any  sensing  system 
should  meet,  and  some  candidate  sensors.  Section  4 deals  with 
the  particular  propagational  and  implementation  features  associated 
with  laser  systems,  and  Section  5 presents  the  system  concept. 

Some  of  the  results  derived  here  were  given  in  a presentation 
in  November  1976  (Kalafus  and  Hallock,  1976)  . 


2.  WIND  SHEAR  DESCRIPTION 


2.1  BOUNDARY  LAYER  SHEARS 

Low-levcl  wind  shears  result  from  the  interaction  of  synoptic 
scale  flow  with  the  surface  of  the  earth.  The  boundary  layer  can 
be  segmented  into  a constant  wind  direction  layer  and  a wind  turn- 
in};  layer.  The  constant  wind  direction  layer  extends  from  the 
surface  up  to  about  150  meters  (500  feet)  and  wind  shear  results 
from  the  boundary  condition  that  the  steady-state  wind  speed  must 
be  zero  at  the  ground.  The  turning  layer  occurs  above  the  con- 
stant wind  direction  layer  and  exhibits  a marked  turning  of  the 
steady-state  wind  vector  with  an  increase  in  altitude. 

Neutral  and  unstable  conditions  rarely  produce  wind  shears 
that  are  hazardous  to  aviation.  Stable  conditions,  on  the  other 
hand,  can  lead  to  appreciable  wind  shears.  Negative  buoyant 
forces  suppress  the  turbulence  and  decouple  the  layers,  allowing 
them  to  slip  relative  to  each  other.  With  very  stable  conditions 
the  flow  is  laminar,  very  little  mixing  occurs,  and  the  overriding 
air  remains  separated  from  the  underlying  air.  Large  changes  in 
wind  direction  (in  excess  of  40®)  and/or  wind  magnitude  (15m/sec, 
or  30  knots)  are  possible  in  altitude  changes  of  less  than  30 
meters  (100  feet).  This  is  discussed  further  in  Appendix  C. 

During  the  night  an  inversion  can  build  up  and  relax  the  con- 
straint imposed  on  the  wind  by  daytime  mixing.  As  a result  the 
wind  at  the  top  of  the  inversion  accelerates,  becomes  supergeo- 
stropic,  and  oscillates  inert ial ly- -a  low-level  nocturnal  jet. 

The  low-level  jet  is  a thin,  well-defined  region  of  high-speed 
air  at  typical  altitude  of  about  300  meters,  with  no  indication  of 
such  a wind  near  the  ground.  The  jet  begins  to  build  up  in  the 
late  afternoon  as  the  earth  cools,  reaches  its  maximum  in  the 
middle  of  the  night,  and  decays  in  the  early  morning.  The  noctur- 
nal jets  can  produce  significantly  large  values  of  low-level  wind 
shear.  At  the  peak  of  the  jet  the  winds  in  its  core,  between  250 
and  700  meters  (800  and  2200  feet)  in  altitude,  can  attain  between 
15  and  35  m/sec  (30  to  70  knots),  decreasing  to  5 or  10  m/sec  (10 


or  20  knots)  between  1000  and  1400  meters  (3000  and  4500  feet) , 
and  to  0 at  the  ground  (see  Fig.  2-1).  It  is  a local  phenomenon 
which  can  be  1500  kilometers  long  (800  nm)  and  70  to  700  kilometers 
wide  (40-400  nm) . In  the  United  States  the  occurrence  of  the 
nocturnal  jet  is  most  pronounced  in  the  Great  Plains.  Tower 
measurements  have  recorded  low-leveTjet  shears  in  excess  of  0.4 
m/sec/m  (24  kphf ) . This  corresponds  to  a headwind  change  of  12 
meters/second  (23  knots)  within  a change  of  altitude  of  30  meters 
(100  feet)  (a  corresponding  change  of  2-3  meters/second  (4-6 
knots)  can  be  hazardous  for  some  aircraft) . 

A vertical  looking  probe  will  adequately  measure  such  a shear. 
These  conditions  occur  over  large  areas  and  last  for  long  periods 
of  time.  Thus  the  probe  measurements  are  highly  correlated  in 
space  and  time  with  aircraft  flight  profiles  so  that  the  measured 
information  should  provide  the  pilot  with  an  accurate  picture  of 
what  conditions  to  expect. 

This  type  of  condition  is  less  dangerous  to  aircraft  than 
fronts  or  storms  for  several  reasons:  (1)  airport  anemometers 
give  good  indications  of  wind  speeds  at  the  surface,  so  that 
pilots  have  some  indication  of  the  presence  of  a shear  by  compari- 
son of  reported  surface  winds  and  winds  aloft;  (2)  reports  of 
previous  pilots  tend  to  be  consistent  and  reliable;  (3)  the  shear 
at  the  upper  part  of  the  jet  is  at  a high  enough  altitude  to  be 
manageable;  (4)  there  are  fewer  operations  during  this  period  of 
time;  and  (5)  visibility  is  good,  and  visual  cues  are  clear. 

2.2  FRONTAL  SYSTKMS 

Fronts  are  the  zones  of  discontinuity  between  cold  and  warm 
air,  typically  polar  and  tropical  air  masses.  These  zones  are 
characterized  by  differences  in  temperature,  wind,  and  humidity 
between  the  adjacent  air  masses,  and  are  frequently  accompanied 
by  stormy  weather.  The  effect  of  the  earth's  surface  is  to  slow 
the  advance  of  a front,  so  that  warm  air  always  lies  above  the 
cold  air  of  the  front.  This  frictional  effect  also  causes  a cold 
front  to  have  a "nose"  so  that  the  slope  of  a rapidly  advancing 


/ U- 

, u;  w 

/ CL  H 

I O 


/Ain.i 


o 

I/) 


133d  ‘3aniinv 


FIGURE  2-1.  VERTICAL  WIND  PROFILE  ASSOCIATED 
WITH  A NOCTURNAL  JET 
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cold  front  can  be  quite  steep,  vertical,  or  even  "bent  over"  in 
the  lowest  300  meters  (1000  feet)  (Brundidge,  1965).  A warm 
front,  on  the  other  hand,  has  a very  shallow  slope  near  the  ground, 
usually  less  than  a degree  (see  Figure  2-2).  On  the  synoptic 
scale,  the  wind  changes  across  a front  are  primarily  that  of  wind 
direction  rather  than  speed.  In  the  Northern  Hemisphere  these 
wind  changes  are  characterized  by  a "left  turn"  at  the  front 
(Figure  2-2).  This  does  not  always  hold  true  on  a smaller  scale, 
because  localized  changes  in  direction  and  speed  may  occur.  The 
zones  of  transition  vary  from  1-300  kilometers  (0.5-160  nm)  in 
width.  When  the  zone  width  is  narrow,  sharp,  and  potentially 
hazardous  changes  in  wind  occur. 

A fast-moving  cold  front  is  usually  accompanied  by  hazardous 
flying  weather,  including  scattered  thunderstorms  when  the  warm 
air  is  moist.  Even  under  dry  conditions,  the  wind  changes  in  the 
frontal  zone  can  be  abrupt  and  large.  Since  an  aircraft  experi- 
ences a passage  through  such  a vertical  zone  as  a sudden  change 
in  headwind  (horizontal  shear)  , the  effect  is  much  the  same  as  for 
an  aircraft  descending  through  an  inversion  (vertical  shear).  It 
can  be  seen  from  Figure  2-2  that  a vertical  - looking  probe  would 
not  detect  the  presence  of  a cold  front  until  it  reached  the  air- 
port. Since  slopes  of  cold  fronts  vary  widely,  the  time  record  of 
such  a probe  would  not  be  able  to  allow  estimates  of  frontal 
speed  or  severity  of  danger  to  an  approaching  or  departing  aircraft. 

Warm  fronts  are  also  associated  with  abrupt  wind  shifts;  the 
Iberia  Air  Lines  accident  at  Logan  (NTSB,  1974c)  appears  to  have 
been  associated  with  a warm  front  wind  shear.  Warm  fronts  have 
shallow  slopes  near  the  ground,  typically  less  than  one  degree, 
so  that  approaching  and  departing  aircraft  operating  at  angles  of 
3-6  degrees  experience  a warm  front  as  a vertical  shear,  much 
like  an  inversion  shear.  A vertical -looking  probe  would  detect 
the  presence  and  measure  the  intensity  of  a warm  front  shear,  but 
would  not  give  speed  or  slope  (see  Fig.  2-3). 
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FIGURE  2-2.  TYPICAL  FRONTAL  PROFILES 


FIGURE  2-3.  WARM  FRONT  DETECTION 
WITH  A VERTICAL  PROBE 


2 . 3 THUNDERSTORMS 

In  mature  thunderstorms  both  updrafts  and  downdrafts  occur 
which  make  it  difficult  for  an  approaching  aircraft  to  hold  the 
glidepath.  In  the  lowest  300  meters  (1000  feet)  of  altitude,  the 
effect  of  the  earth’s  surface  is  to  cause  an  outflow  from  a down- 
draft.  Usually  a downdraft  is  accompanied  by  heavy  rain,  but 
there  have  been  cases  where  the  precipitation  evaporated  before 
reaching  the  ground  (NTSB,  1976b)  . The  outflow  results  in  a gust 
front  ahead  of  the  downdraft  of  a thunderstorm  (Figure  2-4)  • The 
gust  front  may  precede  the  precipitation  by  up  to  18  nm  (30  km) . 
The  gust  front  varies  from  1-3  kilometers  (0.7-2  nm)  in  thickness 
I (Goff,  1975) . An  aircraft  penetrating  the  gust  front  usually 

experiences  updrafts,  an  increasing  headwind,  and  a change  in 
crosswind.  The  slope  of  the  front  is  much  steeper  than  the  path 
of  aircraft,  so  that  it  is  characterized  by  a horizontal  shear. 
Headwind  gradients  reported  by  Goff  (1975)  reached  values  of 
0.1  m/sec/m  (6  knots  per  thousand  feet),  which  corresponds  to  a 
headwind  speed  change  of  3.5  m/sec  (7  knots)  in  5 seconds  for  an 
aircraft  traveling  at  70  m/sec  (140  knots). 

Within  a large  thunderstorm  system  there  can  be  several  in- 
tense downburst  cells  (Caracena,  1976  and  Fujita,  1975).  These 
cells  can  produce  downdrafts  of  10  m/sec  (33  fps)  and  extreme  wind 
speed  changes.  In  the  Stapleton  accident  (NTSB,  1976b)  a loss  of 
headwind  of  12  m/sec  (23  knots)  in  5 seconds  was  recorded.  'The 
cells  are  as  little  as  3 kilometers  (2  nm)  across,  and  can  be 
traversed  by  an  aircraft  in  less  than  a minute.  Similarly,  behind 
the  gust  front,  large  scale  turbulence  can  take  the  form  of  up- 
drafts and  downdrafts  or  secondary  surges;  these  secondary  surges 
exhibit  similar  characteristics  to  the  downburst  cells,  but  are 
not  accompanied  by  rain,  and  are  not  visible  to  the  pilot.  Down- 
draft  cells  such  as  these  are  believed  to  be  the  most  dangerous 
hazard  to  landing  and  departing  aircraft.  Within  40  seconds  the 
pilot  may  have  to  successively  handle  an  increasing  headwind  (gust 
■ front),  decreasing  headwind,  blinding  rain,  downdraft  and  cross- 

wind,  increasing  tailwind,  and  increasing  headwind  (gust  front). 


FIGURE  2-4.  GUST  FRONT  OF  A THUNDERSTORM 
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while  simultaneously  acquiring  visual  contact  with  the  runway 
This  situation  is  portrayed  in  Figure  3-4  and  Table  3-2. 


2-9/2-10 


r 


I 


I 

i 

i 


L 


'«r 


3.  MAJOR  REQUIREMENTS  FOR  A GROUND-BASED  SYSTEM 

Within  the  constraints  of  cost,  system  requirements  are 
dictated  primarily  by  the* intended  use  of  the  information  acquired. 
In  the  case  of  wind  shear,  the  data  will  be  used  to  detect  the 
existence  of  a hazardous  shear  event,  and  provide  the  pilot  and 
controller  with  sufficient  information  on  the  intensity  and  alti- 
tude of  the  expected  encounter,  so  that  decisions  can  be  safely 
made  relative  to  (1)  execution  of  a missed  approach  (or  choosing 
not  to  depart)  or  (2)  continuation  of  the  operation  with  a knowl- 
edge of  what  preventive  and  corrective  actions  will  ensure  the 
safety  of  the  flight.  In  the  present  context,  this  means  detecting 
shear  zones,  tracking  them,  measuring  their  intensity  and  estimat- 
ing their  time  of  arrival  at  the  corridor. 

3.1  SYSTEM  RANGE 

System  range  should  be  large  enough  to  protect  the  approach 
and  departure  corridors  out  to  a point  where  an  aircraft  has  a high 
probability  of  recovering  after  encountering  a severe  wind  shear 
event.  The  National  Transportation  Safety  Board  has  made  the 
following  recommendations  (NTSB,  1976b): 

"Expedite  the  program  to  develop  and  install  equipment  which 
would  facilitate  the  detectior  and  classification,  by  severity,  of 
thunderstorms  within  4 nm  of  the  departure  or  threshold  ends  of 
active  runways  at  airports  having  precision  instrument  approaches. 
(A-76-32.)" 

"Install  equipment  capable  of  detecting  variations  in  the 
speed  of  the  longitudinal,  lateral,  and  vertical  components  of  the 
winds  as  they  exist  along  the  projected  takeoff  and  approach  flight- 
paths  within  1 nm  of  the  ends  of  active  runways  which  serve  air 
carrier  aircraft.  (A-76-33.)" 

On  the  basis  of  these  and  other  recommendations,  the  FAA  Wind 
Shear  Program  Office  has  adopted  the  requirement  of  coverage  to  the 
outer  marker,  until  results  from  simulations,  hazard  and  accident 
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studies  indicate  otherwise.  This  implies  a system  range  of  about 
7 nm  (13  km)  for  a sensor  placed  in  a centerfield  location. 

Rather  than  deriving  a specific  number  for  system  range,  the 
remainder  of  this  discussion  is  aimed  at  outlining  the  steps  by 
which  this  number  can  be  determined,  using  results  from  the 
various  study  efforts.  The  major  considerations  are: 

a.  Where  along  the  flight  path  have  shear  situations  been 
encountered  which  historically  have  caused  (or  contri- 
buted to)  accidents  (or  near  accidents)? 

b.  Based  on  simulations  and  reasonable  operating  assumptions, 
what  is  the  minimum  coverage  such  that,  if  no  shear  events 
exist  in  the  coverage  area,  the  probability  of  an 
accident  is  small? 

c.  How  do  these  factors  relate  to  the  type  of  shear  en- 
countered? 

The  accident  analysis  should  be  used  to  back  up  the  coverage 
requirement  resulting  from  (b)  and  (c) . 

The  first  question  to  be  addressed  is  whether  system  range 
should  be  based  on  approach  requirements  or  departure  requirements. 
On  takeoff,  the  aircraft  is  typically  climbing  at  rates  of  1200 
fpm  (6  m/sec)  or  more,  and  reaches  a safe  airspeed  and  altitude 
within  0.5-1  nm  (1-2  km)  of  the  runway  end.  On  approach,^  the 
aircraft  is  operating  near  stall  speeds  farther  out  from  touch- 
down. Thus  system  range  must  be  premised  on  approach  requirements. 
System  range  is  determined  by:  first,  establishing  the  minimum 
range  for  which  an  aircraft  can  readily  recover  after  encountering 
an  event  farther  out  (critical  zone);  second,  establishing  th 
amount  of  advance  warning  time  needed  to  assure  that  no  sheai 
event  can  move  into  the  critical  zone  before  pilots  can  be 
alerted  . 

Clearly,  a shear  which  causes  a plane  to  stall  is  dangerous 
at  any  altitude.  There  are  no  accidents  of  this  type  known  to  the 
author.  However,  shears  can  be  dangerous  without  causing  stalls; 
the  extent  of  the  danger  then  depends  on  the  altitude  of  the 
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aircraft  at  the  time  of  encounter- -the  lower  the  altitude  the  more 
dangerous  the  situation.  The  zone  to  be  protected  is  the  approach 
corridor  from  threshold  out  to  a distance  sufficiently  far  so  that, 
if  an  aircraft  encounters  a shear  event  outside  the  zone, 
recovery  is  virtually  certain  (see  Figure  3-1). 

The  ’bility  of  an  aircraft  to  recover  depends  on: 

a.  The  mass,  available  thrust,  and  responsiveness  of  the 
aircraft . 

b.  The  type  of  shear  event,  its  intensity,  and  its  duration. 

c.  The  information  available  to  the  pilot. 

d.  The  pilot's  skill. 

A ground-based  system  should  be  conservative:  the  aircraft 
studied  should  be  primarily  high-mass,  slow-response  aircraft; 
the  shear  events  assumed  should  be  based  on  the  worst  observed 
conditions:  minimal  information  should  be  assumed  (no  heads-up 

display  or  inertial  system,  for  examples),  and  pilot  response 
should  be  assumed  to  be  slow. 

There  are  several  reasons  for  believing  that  the  danger 
decreases  for  events  encountered  farther  up  the  glide  path: 

a.  The  aircraft  has  more  altitude  and  time  for  recovery; 

b.  The  pilot  has  fewer  duties:  decision  height  is  not  yet 
reached,  the  aircraft  may  still  be  on  autopilot  and 
margins  are  greater. 

c.  In  the  important  case  of  a downburst  event,  the  danger 
is  greatest  below  500  feet;  this  is  discussed  at  some 
length  below. 

Collectively,  these  factors  imply  a marked  decrease  in  the 
danger  from  a shear  event  encountered  above  about  500-feet  alti- 
tude (ISO  m) . 

The  downdrafts  associated  with  a thunderstorm  constitute  the 
single  most  dangerous  wind  shear  event  known.  Convincing  evidence 
has  only  recently  been  obtained  on  the  severity  of  these  events, 
but  this  evidence  is  dramatic. 
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THE  PROTECTED  'CRITICAL  ZONE' 
SHOULD  BE  LARGE  ENOUGH  THAT  IF 
AN  AIRCRAFT  ENCOUNTERS  A SEVERE 
SHEAR  EVENT  BEFORE  REACHING  IT, 
RECOVERY  IS  HIGHLY  PROBABLE. 


FIGURE  3-1,  APPROACH  CORRIDOR 
REQUIRING  PROTECTION 
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While  downdrafts  from  a thunderstorm  outflow  can  occur  over 
a larger  area,  such  wide-scale  drafts  are  believed  to  be  small, 
on  the  order  of  2-5  feet  per  second  (0.6-1. 5 m/sec)  at  altitudes 
below  500  feet  (150  m) . In  downburst  cells  the  intense  drafts  are 
tightly  confined,  being  on  the  order  of  1000-2000  feet  (300-600  m) 
wide.  This  statement  is  premised  on  limited  data,  but  the  effect 
is  observed  in  data  from  Goff  (1976),  Caracena  (1976)  and  Fujita 
(1975).  Below  300-500  feet  (90-150m)  the  effect  of  the  ground 
is  to  convert  the  downwash  into  an  outflow  (see  Figure  3-2).  The 
outflow  is  approximately  radially  outward  from  the  center  of  the 
downwash  cell;  the  magnitude  varies  substantially,  as  the  wind 
apparently  tends  to  concentrate  in  sectors.  Figure  3-3,  taken 
from  Caracena  (1975),  shows  the  horizontal  streamlines  associated 
with  a thunderstorm  outflow.  The  isotachs  in  that  study  indicate 
wind  speeds  that  vary  from  6 knots  to  40  knots  in  the  radial 
streamlines  around  the  downwash  cell  associated  with  the  Denver 
Stapleton  accident  (NTSB-,  1976b)  . An  aircraft  traversing  the  cell 
below  300-500  feet  (90-150  m)  would  encounter  a decreasing  head- 
wind shear,  coupled  with  a downdraft  and  possible  crosswinds.  The 
archetypal  model  is  shown  in  Figure  3-4  and  Table  3-2. 

To  ascertain  the  relative  effects  of  the  downdraft  and  a 
horizontal  shear  associated  with  a downburst  cell,  the  analog 
computer  study  by  Snyder  (1968)  was  used.  There  he  simulated  a 
swept-wing  transport  aircraft  on  final  approach,  and  subjected  it 
to  several  events:  (1)  a sudden  5 knot  (500  fpm,  or  2.5  m/sec) 
downdraft,  sustained;  (2)  a sudden  15  knot  (8  m/sec)  drop  in  air- 
speed, sustained;  and  (3)  a sudden  shear  of  5 (and  10)  kphf 
(0.084  and  0.17  m/sec/m),  also  sustained.  Controls  are  fixed, 
providing  a good  picture  of  what  happens  in  the  absence  of  pilot 
intervention.  The  response  to  a downdraft  is  a temporary  pitchup, 
which  increases  lift,  thus  countering  the  downward  thrust  of  the 
draft  to  some  extent;  eventually  the  aircraft  stabilizes  in  the 
new  air  mass  at  a sink  rate  equal  to  the  original  sink  rate  plus 
the  downdraft  speed.  The  response  to  a headwind  loss  is  an 
initial  pitchdown;  this  and  the  loss  of  airspeed  both  act  to 
increase  the  sink  rate.  Snyder's  computations  were  premised  on 
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a rather  fast  approach  speed  of  160  knots  (95  m/sec  airspeed)  and 
initial  trim  at  an  800  fpm  (4  m/sec)  sink  rate.  Table  3-1  shows 
the  sink  rates  associated  with  the  three  shear  types. 


TABLE  3-1.  SINK  RATES  ASSOCIATED  WITH  SHEAR 
EVENTS  (FROM  SNYDER,  1968) 


SHEAR  EVENT 

AFTER  5 SECONDS 

AFTER  10  SECONDS 

5 knot  downdraft 

1000  fpm  (5  m/sec) 

1300  fpm  (6.5  m/sec) 

15  knot  headwind  loss 

1900  fpm  (9.5m/sec) 

2500  fpm  (12.5m/sec) 

5 kphf  shear 

950  fpm  (4.8m/sec)  . 

1400  fpm  (7  m/sec) 

10  kphf  shear 

1100  fpm  (5.5m/sec) 

2100  fpm  (10.5m/sec) 

Using  Snyder's  formulation,  TS,C  performed  computations  which 
duplicated  his  results  and  were  run  for  other  shear  severities 
and  approach  speeds  of  145  and  120  knots.  The  results  were  con- 
sistent with  Snyder's. 

It  can  be  seen  from  the  table  that  a 5 knot  downdraft  is 
comparable  to  a 5 kphf  (0.08  m/sec/m)  shear  for  about  10  seconds; 
thereafter  the  shear  is  worse,  since  the  sink  rate  continues  to 
increase.  Whether  a 10  knot  downdraft  is  comparable  to  a 10  kphf 
(0.17  m/sec/m)  shear  is  not  known,  but  it  appears  to  be  a reason- 
able assumption.  There  is  no  firm  relationship  between  the  down- 
draft  speed  and  the  equivalent  vertical  shear  in  a downburst  cell, 
for  the  obvious  reason  that  the  air  mass  flow  is  not  uniformly 
distributed  in  angle.  Thus  the  shear  encountered  by  an  aircraft 
depends  on  the  direction  from  which  the  cell  is  approached,  as 
well  as  the  region  within  the  cell  that  is  traversed.  However, 
based  on  a simple  flow  model,  it  appears  that  a 5 kphf  (0.08  m/ 
sec/m)  shear  is  a typical  average  shear  that  would  be  associated 
with  a 5 knot  (2.5  m/sec)  downdraft. 

The  conclusion  to  be  drawn  from  this  is  that  in  general  an 
aircraft  encountering  a downdraft  at  higher  altitudes  will  deviate 
from  the  glidepath  less  than  if  it  encounters  the  same  downburst 
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FIGURE  3-2.  DOWNBURST  CELL  SHOWING  OUTWASH 
DUE  TO  THE  PRESENCE  OF  THE  GROUND 
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FIGURE  3-3.  STREAMLINES  ASSOCIATED  WITH 
STAPLETON  ACCIDENT,  AUGUST  197S 
(FROM  CARACENA,  1976) 


CLASSIC  DOWNBURST  ENCOUNTER- -STREAMLINES 
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LIGHT  HEADWIND  CLIMBS  OUT  SWEAT 


at  lower  altitudes  where  the  downdraft  is  coupled  with  a diminish- 
ing headwind  shear. 

The  extent  of  the  critical  zone  does  not,  by  itself,  estab- 
lish the  system  range.  It  was  asserted  earlier  that  once  the 
extent  of  the  critical  zone  is  established,  the  system  should 
provide  advance  warning.  Thus  the  system  coverage  must  be  large 
enough  to  guarantee  that  an  event  does  not  move  into  the  critical 
zone  before  it  can  be  detected,  identified,  and  a warning  issued 
and  acknowledged.  This  is  discussed  in  the  next  section. 

3.2  DATA  UPDATE  RATE  AND  COVERAGE 

The  maximum  rate  at  which  the  airport  wind  shear  situation 
is  surveyed  is  limited  by  the  sensor.  Regardless  of  the  sensor 
used,  the  minimum  update  rate  depends  on  the  coverage,  the  dura- 
tion and  speed  of  the  phenonienon,  the  time  required  to  divert  an 
aircraft  from  landing  or  departing,  and  the  basic  concept  of  the 
use  of  the  information.  Since  there  is  frequently  a cost/perform- 
ance tradeoff  between  update  rate  and  other  system  parameters, 
the  minimum  update  rate  is  the  important  system  parameter  to  be 
determined. 

Minimum  coverage  and  update  rate  are  determined  by  landing, 
rather  than  departure  requirements,  because  the  approaching  air- 
craft is  slower,  and  is  in  a more  critical  conf iguratipn  near  the 
ground  for  a longer  period  of  time.  If  the  meteorological  situa- 
tion did  not  change  during  the  approach,  it  would  be  sufficient 
to  protect  the  zone  shown  in  Figure  3-1.  That  is,  by  assuming 
that  there  was  no  shear  event  in  the  critical  zone,  approaching 
aircraft  would  be  protected. 

This  kind  of  protection  would  be  sufficient  for  wide-scale 
synoptic  events  such  as  those  encountered  in  inversions  and  shallow 
sloped  warm  fronts.  There  the  basic  wind  patterns  change  very 
slowly  with  time.  Where  thunderstorms  or  cold  fronts  with  steep 
slopes  are  involved,  however,  wind  patterns  can  change  significant- 
ly in  less  than  a minute.  In  these  situations  the  information 
must  be  updated  frequently,  and  the  coverage  must  be  expanded  to 
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detect  shear  events  moving  toward  the  corridor. 

An  examination  of  gust  front  speeds  (Goff,  1975  and  1976) 
and  storm  system  motions  (Caracena,  1976  and  Fujita,  1975) 
suggests  that  velocities  of  10-15  knots  are  typical,  while 
velocities  of  30  knots  have  been  observed.  There  were  two  cases 
in  Goff  (1975)  of  faster  gust  front  motion,  but  the  associated 
shears  were  mild.  For  the  purposes  of  this  analysis,  30  knots 
(15  m/sec)  will  be  used  as  a working  figure  for  rapidly  moving 
shear  events.  In  the  2-2.5  minutes  that  an  approaching  aircraft 
takes  to  fly  from  the  outer  marker  to  touchdown,  such  a wind 
shear  event  can  move  over  a nautical  mile. 

Any  real  sensor  has  a time  delay  associated  with  it:  the 
effect  must  be  observed  for  some  period  of  time  before  it  registers; 
a sampling  system  has  a delay  between  zero  and  the  update  period, 
when  it  first  appears  to  the  sensor.  Thus  any  point  sensor  or 
vertical-looking  sensor  will  not  register  an  alarm  until  the 
wind  shear  event  has  been  present  above  the  sensor  for  some  period 
of  time.  If  this  period  of  time  is  excessive,  e.g.,  more  than  a 
minute,  such  a sensor  located  in  the  critical  zone  of  Figure  3-5 
would  not  provide  adequate  protection.  Even  half  a minute  would 
be  marginal,  since  the  shear  event  could  move  about  1500  feet 
(450  m) ; shear  events  can  occur  on  a scale  of  1000-1500  feet 
(300-450  m)  (see  section  3.3)  and  thus  could  be  outside  the 
corridor  at  one  time  and  inside  the  corridor  half  a minute  later. 

Where  sensor  coverage  encompasses  a wider  area,  the  shear 
event  can  be  detected  before  it  reaches  the  critical  zone.  If 
this  area  is  sampled  sufficiently  often,  the  event  can  be 
observed  several  times  before  it  reaches  the  operating  corridor. 

It  is  anticipated  that  a downburst  cell  will  maintain  its  shape 
and  intensity  to  a sufficient  degree  to  be  tracked;  i.e.,  be 
unique  and  consistent  enough  from  sample  to  sample  so  that  its 
intensity,  speed,  and  direction  can  be  determined.  This  would 
enable  an  algorithm  to  be  employed  which  would  give  advance 
warning  of  the  time  of  arrival  of  the  event  in  the  operating 
corridor. 
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A hit  is  defined  as  one  observation  of  a shear  event.  Two 
hits  are  necessary  to  establish  a velocity,  and  probably  3 to  4 
hits  will  prove  necessary  for  validation  and  elimination  of  false 
alarms,  before  an  alarm  is  issued.  The  additional  coverage  this 
implies  is  shown  in  Figure  3-5  for  3-  and  4-hit  criteria,  assuming 
a fast  30  knot  front  velocity.  Also  included  is  the  time  T 
required  between  the  time  an  alarm  is  issued  and  the  time  the 
pilot  is  alerted,  or  a decision  made  to  close  the  runway. 

The  establishment  of  the  time  duration  T is  an  important 
parameter  which  needs  to  be  addressed.  Once  it  is  established, 
the  additional  coverage  required  for  each  approach  corridor  would 
be  that  shown  in  Figure  3-6.  The  composite  coverage  for  a typical 
airport  would  then  resemble  Figure  3-7. 

3.3  SPATIAL  RESOLUTION 

Spatial  resolution  depends  on  the  size  of  the  phenomenon: 
the  thickness  of  the  shear  zone,  its  shape,  and  its  extent.  It 
also  depends  on  the  aircraft  response  scale.  Even  if  shears 
existed  on  a 10  foot  (3  m)  scale,  if  an  aircraft  responded  on  a 
100  foot  (30  m)  scale  it  would  not  be  affected  by  them;  stated 
another  way,  shear  effects  which  are  seen  by  an  aircraft  at  high 
frequencies  compared  with  the  frequency  response  of  the  aircraft 
are  not  significant. 

A typical  aircraft  will  deviate  from  the  flight  path  in 
response  to  frequencies  below  about  one  radian  per  second. 
Translated  to  an  aircraft  traveling  120  knots  (62  m/sec),  this 
corresponds  to  a horizontal  scale  of  about  1300  feet  (400  m) , 
and  a vertical  scale  of  70  feet  (20  m) . This  suggests  that  a 
horizontal  resolution  of  1500  feet  (450  m)  or  less  would  be 
desirable  from  the  point  of  view  of  aircraft  response.  Similarly, 
a vertical  resolution  of  75  feet  (23  m)  or  less  would  be  desirable. 

This  discussion  is  included  to  show  that  the  relative  scales 
are  reasonable:  disturbances  on  a smaller  scale  than  1000  feet 
(300  m)  may  cause  aircraft  control  surface  motion  and  buffeting, 
but  do  not  move  the  aircraft  significantly.  Of  equal  importance 
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is  the  scale  on  which  wind  changes  take  place  in  shear  events. 

If  these  changes  occur  only  on  scales  much  larger  than  1000  feet, 
resolution  would  not  need  to  be  so  fine-grained. 

For  large-scale  synoptic  phenomena,  wind  gradients. are 
gradual  in  the  horizontal  dimension,  and  are  significant  only 
in  the  vertical  dimension.  Thus,  the  shear  events  are  adequately 
represented  by  a vertical  shear  measurement  such  as  would  be 
obtained  by  a vertical-looking  probe. 

In  the  case  of  the  boundary  layer  shears,  a logarithmic 
profile  near  the  ground  yields  a good  approximation  to  the 
observed  values.  Here  the  winds  change  in  a regular  fashion 
between,  say,  500  feet  (150  m)  height  AGL  and  the  surface.  A 
height  resolution  of  100  feet  (30  m)  would  be  adequate  to  describe 
the  shear  conditions. 

However,  the  profiles  found  in  fronts  (NTSB,  1974c  and  Coons, 
1976)  often  exhibit  less  regular  profiles,  ones  whose  winds 
change  significantly  over  a 100  foot  (30  m)  altitude  movement, 
but  change  less  sharply  elsewhere.  In  these  cases  a 100  foot 
resolution  would  not  be  fine  enough  to  characterize  the  shear; 
this  is  shown  in  Figure  3-8.  On  the  other  hand,  a 50  foot  (15  m) 
resolution  is  quite  adequate.  Thus  for  vertical -looking  probes, 
a 50  foot  altitude  resolution  is  recommended. 

Cold  fronts  and  thunderstorm  downwash  surges  are  not  as  well 
characterized  by  vertical  shear  derived  from  a vertical  profile 
measurement  as  by  horizontal  shear.  An  aircraft  on  final  approach 
experiences  both  vertical  and  horizontal  shears  as  changes  in 
headwind  (actually,  in  crosswind  as  well),  so  there  are  similarities 
between  the  two  types.  However,  with  different  descent  rates 
aircraft  flying  through  a horizontal  shear  will  experience  similar 
headwind  changes. 

Cold  fronts  have  "noses"  caused  by  frictional  retardation 
at  the  ground.  The  slope  near  the  ground  can  be  quite  shallow 
(slow-moving),  quite  steep,  or  even  "bent  over."  The  scale  is 
such  that  the  shear  could  be  horizontal  or  vertical,  or  something 
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in  between.  Brundidge  (1965)  reports  examples  where  the  wind 
shift  zones  vary  from  1000  to  5000  feet  in  extent  (300  m to 
1500  m) . 

Gust  fronts  associated  with  thunderstorm  outflows  typically 
exhibit  updrafts  and  severe  shears  which  an  aircraft  would 
experience  as  an  increasing  headwind  as  it  travels  through  them. 
An  example  is  shown  in  Figure  3-9,  taken  from  storm  data  at  NSSL  . 
(Goff,  1975).  The  contours  in  the  bottom  record  represent 
constant  values  of  the  wind  component  normal  to  the  front.  There 
is  a shear  zone  about  a mile  long  wherein  the  headwind  gradient 
experienced  by  an  aircraft  having  a ground  speed  of  120  knots 
(62  m/sec)  would  be  about  16  knots  (8  m/sec)  per  10  seconds  of 
travel,  and  a total  headwind  increase  of  32  knots  (16  m/sec)  in 
about  45  seconds.  An  interesting  aspect  of  this  is  that  an 
aircraft  will  experience  an  increasing  headwind  regardless  of 
which  way  it  flies  through  the  gust  front.  For  wind  speed  con- 
tours like  this  one,  even  range  resolution  as  large  as  4000  feet 
would  satisfactorily  characterize  the  total  wind  change,  since 
on  either  side  of  the  gust  front  fairly  uniform  winds  are  en- 
countered. This  large  a resolution  would  not  adequately  estimate 
the  shear  for  reasons  described  in  Figure  3-8,  but  1000-2000 
foot  (300-600  m)  resolution  is  quite  satisfactory. 

Intense,  concentrated  downwash  cells  are  believed  to  be 
associated  with  the  Eastern  66  and  Continental  426  accidents 
(Fujita,  Caracena) . These  cells  appear  to  be  as  little  as  2nm 
(4  km)  in  diameter,  and  are  accompanied  by  severe  downdrafts. 

The  headwinds  encountered  by  aircraft  traversing  the  cell  would 
resemble  those  depicted  in  Figure  3-10.  Along  paths  AA  and  A'A' 
cross-track  gusts  would  also  be  experienced,  complicating  pilot 
workload  and  causing  lateral  deviations  from  the  extended  runway 
centerline.  Downdrafts  would  be  experienced  in  all  paths,  being 
worst  in  path  BB.  Headwind  changes  observed  in  such  cells  have 
ranged  from  20  knots  (10  m/sec)  in  12  seconds  (Coons  and  Mandel, 
1976)  to  40  knots  (20  m/sec)  in  5 seconds  (NTSB,  1976b);  in  the 
first  case,  the  headwinds  changed  from  zero  to  25  knots  (13  m/sec ] 


3-19 


b.  PLAN  VIEW 


HEADWIND  t 


c.  WIND  PROFILE 
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to  zero  in  300  feet  (90  m)  of  altitude  change,  i.e.,  in  about  25 
seconds.  This  translates  to  about  one  mile  of  travel.  For 
reasons  already  stated,  resolution  sizes  larger  than  about  2000 
feet  (600  m)  would  not  adequately  characterize  the  wind  shear, 
and  1000  feet  (300  m)  would  be  preferable. 

Thus,  based  on  the  meteorological  phenomena  and  on  air- 
craft response  times,  1000  feet  (300  ra)  is  the  preferred  resolu- 
tion, while  1500  feet  (450  m)  is  probably  acceptable.  If  finer- 
grained  resolution  could  be  obtained  without  sacrificing  other 
system  parameters,  it  would  be  desirable,  but  not  necessary. 

The  cross-track  resolution  does  not  appear  to  be  critical, 
but  it  should  be  fine  enough  to  enable  a reasonably  accurate 
assessment  of  the  width  of  the  shear  zone,  and  the  variation  of 
intensity  within  it.  From  this  point  of  view,  a resolution  of 
1000-2000  feet  (300-600  m)  is  probably  adequate. 

3 . 4 ACCURACY 

The  wind  shear  accuracy  is  premised  on  the  requirement  of 
achieving  an  acceptably  low  occurrence  of  "false  alarms"  and 
"misses."  A "false  alarm"  occurs  if  the  system  declares  the 
existence  of  a shear  event  when  in  fact  there  is  none.  If  the 
number  of  false  alarms  is  excessive,  the  system  would  not  be 
relied  upon  by  pilots  and  controllers.  A "miss"  is  the  failure 
of  the  system  to  detect  a hazardous  shear--the  impact  of  this 
is  obvious. 

The  issue  of  false  alarms  and  misses  is  more  complex  if 
the  system  attempts  to  assess  the  relative  intensity  of  the  shear, 
e.g.,  utilizing  categories  such  as  Light-Moderate-Strong-Severe. 
(See  Table  3-3.) 


TABLE  3*3.  WIND  SHEAR  DESCRIPTORS,  SIMILAR  TO  ICAO 
RECOMMENDATIONS  (ICAO,  1976) 


DESCRIPTION 

SHEAR  MAGNITUDE 

Light 

0-4  kphf  (0-0.07  m/sec/m) 

Moderate 

5-8  kphf  (0.08  - 0.14  m/sec/m) 

Strong 

9-12  kphf  (0.15  - 0.20  m/sec/m) 

Severe 

13  or  more  kphf  (>0.22  m/sec/m) 

An  example  of  a false  alarm  here  would  be  the  declaration 
of  a strong  shear,  when  in  fact  the  shear  is  moderate.  A miss 
would  be  an  under-assessment  of  the  intensity,  e.g.,  the  declara- 
tion of  a moderate  shear  when  the  shear  is  actually  strong. 

For  the  purpose  of  establishing  the  accuracy  requirement, 
a shear/no  shear  model  will  be  employed.  Appendix  D develops 
the  probabilities  of  misses  and  false  alarms  as  a function  of 
wind  shear  error,  taking  into  account  the  fact  that  severe  shears 
are  less  probable  than  light  shears.  Two  system  thresholds  are 
chosen  for  consideration:  4.5  kphf  (0.076  m/sec/m)  and  7.5 
kphf  (0 . 13m/sec/m)  (for  4.5  kphf,  0-4  kphf  is  considered  "no 
shear,"  and  5 or  more  is  considered  a "shear").  The  results  are 
plotted  in  Figure  3-11.  It  can  be  seen  that  if  the  threshold  is 
set  too  high,  the  miss  probability  increases,  while  the  false 
alarm  probability  drops. 

The  chart  should  be  read  in  the  following  way:  for  example, 
given  a measurement  error  of  a kphf,  the  probability  that  a 
shear  of  5 kphf  or  more  will  be  mistakenly  interpreted  as  a 
no-shear  condition  is  given  by  the  miss  curve  for  a ■ 4.5.  A 
different  question  is  the  following:  Given  that  a shear  is  a 
particular  value,  e.g.,  9 kphf,  what  is  the  probability  of  this  l 
being  mistakenly  interpreted  as  a no-shear  condition?  This  is 
given  in  Figure  3-13.  From  this  it  can  be  seen  that  most  of  the 
misses  of  Figure  3-12  occur  for  shears  just  above  threshold. 
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PROBABILITY 


FIGURE  3-11.  VARIATION  OF  MISS  AND  FALSE  ALARM 
PROBABILITIES  WITH  SYSTEM  THRESHOLD 


MISS  PROBABILITIES 


From  these  figures  it  is  evident  that  to  reliably  identify 
shears  of  6 kphf  (0.10  m/sec/m)  a standard  deviation  of  1-2 
kphf  (0.5-1  m/sec)  would  be  required;  this  gives  a 95%  probabil- 
ity of  correct  detection  at  6 kphf.  Even  so,  a 5 kphf  (0.083 
m/sec/m)  shear  would  have  a high  miss  rate.  This  can,  of  course, 
be  improved  if  the  threshold  is  lowered.  However,  if  the  threshold 
is  lowered,  the  boundary  layer  shear  that  routinely  exists  would 
frequently  be  interpreted  as  a shear  event.  This  common  shear  is 
typically  2-3  kphf  near  the  ground,  but  averages  as  much  as  4 kphf 
in  some  areas  (Nancoo,  1974).  These  values  were  obtained  at  obser- 
vation stations  below  150  feet  AGL;  the  shears  between  100  and  300 
feet  will  be  about  half  these  published  values,  assuming  a 
logarithmic  profile.  Thus  a threshold  value  of  3-5  kphf  is  a 
reasonable  one,  which  implies  that  a standard  deviation  require- 
ment of  1-2  kphf  accuracy  is  necessary  to  reliably  assess  shears 
of  5-7  kphf,  values  which  have  been  observed  in  wind  shear 
accidents  (see  Table  3-4). 
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TABLE  3-4.  ESTIMATED  SHEARS  FOR  SEVERAL 
ACCIDENTS  OR  NEAR- ACCIDENTS 


FLIGHT 

LOCATION 

APPROACH 

SPEED, 

KNOTS 

MAXIMUM 

VERTICAL 

SHEAR, 

KPHF 

MAXIMUM 

HEADWIND 

GRADIENTS, 

KNOTS/SEC 

TOTAL 

HEADWIND 

CHANGE, 

KNOTS/SEC 

EA66 

JFK  NY 

140 

-14* 

-2.8 

-18/8 

EA902 

JFK  NY 

140 

-12* 

-2.0 

-14/8 

C0426 

DEN  CO 

160** 

N.A. 

-8.8 

-40/5 

IA933 

BOS  MA 

150 

+7 

+1.4* 

+21/20 

DL516 

CHA  TN 

130 

5* 

-1.0 

-15/15 

* Equivalent  - see  text 
**  Maximum  takeoff  speed 
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The  preceding  discussion  has  used  kphf,  a measure  of  verti- 
cal shear,  for  convenience;  however,  the  reasoning  applies  as 
well  to  horizontal  shear.  The  equivalent  numerical  value  for  a 
horizontal  shear  event  is  found  hy  considering  the  time  rate  of 
headwind  change  experienced  by  an  aircraft  flying  through  verti- 
cal shear  zone  along  the  glidepath:  a significant  vertical  shear 
of  4.5  kphf  would  be  experienced  as  a headwind  change  of  2.4 
knots  (1.2  m/sec)  in  1000  feet  (300m)  of  travel  along  a 3°  glide 
slope.  An  aircraft  traveling  at  160  knots  (82  m/sec)  would  experi- 
ence a headwind  change  of  0.65  knot  (0.33  m/sec)  per  second:  at 
120  knots  (62  m/sec),  the  headwind  change  would  be  0.5  knot  (0.25 
m/sec)  per  second.  Thus  0.5  knot  per  second  or  more  is  considered 
a significant  headwind  gradient,  and  2.5  knots  (1.3  m/sec)  per 
1000  feet  (300  m)  is  considered  a significant  horizontal  shear. 

The  severity  of  the  shear  event  is  not  only  related  to  the 
shear  intensity,  but  also  to  the  duration  of  the  event  and  the  total 
change  in  wind  speed.  Typical  profiles,  somewhat  oversimplified, 
are  shown  in  Figure  3-13;  the  first  two  profiles  are  the  kind 
that  would  be  encountered  by  an  aircraft  flying  down  through  a 
warm  front,  or  a thunderstorm  initial  gust  front.  Downburst  cell 
profiles  would  be  more  like  that  of  Figure  3-13c. 

The  total  wind  change  is  important  because  the  time  required 
for  the  stabilization  of  an  aircraft  in  the  new  air  mass  field 
d'^oeiids  on  it.  Accuracy  for  the  measure  of  total  wind  change  is 
t‘.  ■ critical  than  it  is  for  the  measurement  of  shear  intensity. 

li  fact  it  is  evident  that  for  profiles  a and  b of  Figure 
3-13  there  exist  some  equivalences  between  shear  intensities  and 
total  changes,  at  least  for  a given  aircraft:  for  example,  a 
shear  of  8 kphf  (0.13  m/sec/m)  with  a total  wind  change  of  20 
knots  (10  m/sec)  might  have  a similar  effect  on  an  aircraft  as  a 
shear  of  12  kphf  (0.20  m/sec/m)  with  a total  wind  change  of  10 
knots  (5  m/sec).  More  work  needs  to  bo  done  in  this  area  because 
it  may  relax  the  accuracy  requirements  and  improve  the  effectiveness 
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b.  INCREASING  HEADWIND  | 


c.  DOWNBURST  CELLS 


FIGURE  3-13.  HEADWIND  PROFILES  ENCOUNTERED  BY  A 
MOVING  AIRCRAFT  FOR  VARIOUS  SHEAR  CONDITIONS 


of  a particular  sensor.  For  the  downburst  profile  of  Figure 
3-13c,  the  situation  is  more  complex,  but  the  same  reasoning 
applies . 


It  should  also  be  pointed  out  that  while  a shear  measure- 
ment accuracy  of  1-2  kphf  (la)  is  a reasonable  svstem  require- 
ment, it  does  not  mean  that  each  sample  of  measurement  must  meet 
this  requirement,  but  that  the  estimate  made  by  a sensing  system 
based  on  the  measurements  should  achieve  this  accuracy. 

3.5  CANDIDATE  SENSORS 

There  are  a number  of  different  ground-based  techniques 
which  have  been  advanced  for  detecting  hazardous  wind  shears. 

Some  of  these  are  discussed  below,  in  terms  of  each  technique, 
its  advantages,  limitations,  and  problems  outstanding.  By 
"technique"  is  meant  the  measurement  method,  data  derived  from 
it,  and  use  of  the  information;  "advantages"  refer  to  the  reasons 
for  using  the  sensor  and  some  commendable  features  of  it;  the 
"limitations"  referred  to  are  inherent  features  of  the  sensor  that 
preclude  it  from  detecting  and  identifying  all  shear  conditions,  as 
well  as  problems  unique  to  the  sensor;  and  by  "problems  outstand- 
ing" are  meant  problem  areas  which  must  be  resolved  if  the  sensor  is 
to  achieve  its  potential. 

Each  technique  is  described,  in  terms  of  its  capabilities, 
etc.,  as  a stand-alone  sensor  system.  Limitations,  in  this  case, 
are  not  meant  as  criticisms  of  the  techniques,  but  merely  to 
delimit  the  expectations  on  its  usefulness.  Use  of  multiple 
sensors,  and  combinations  of  ground  and  airborne  techniques  can 
be  employed  which,  to  varying  degrees  of  success,  meet  the  over- 
all system  requirements. 

a . Anemometers  Near  Runway  Thresholds 

Technique:  Anemometers  are  placed  beside  active  runways  at 
each  end.  Readings  are  relayed  via  radio  link 
or  ground  lines  to  the  tower.  Pilots  are  alerted 
as  to  ground  wind  speed  and  direction  at  the 
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appropriate  runway,  either  by  the  controller  or 
the  ATIS  broadcast.  Departing  aircraft  would  be 
advised  if  the  vector  difference  between  wind 
velocities  at  opposite  ends  of  the  runway  ex- 
ceeded some  threshold  like  15  knots. 

Advantages:  Inexpensive  and  utilizes  familiar  and  relatively 
simple  equipment. 

Limitations : Does  not  detect  shear  events  beyond  airport 

boundaries;  does  not  provide  advance  warning  of 
shear  events  moving  into  the  corrdor;  does  not 
detect  synoptic  shears  above  the  airport. 

Problems  Outstanding:  Establish  whether  tradeoff  between 
false  alarms  caused  by  gusts  and  time  lag 
associated  with  smoothing  is  feasible. 

b.  Anemometers  Beyond  the  Middle  Markers 

Technique:  Anemometers  are  placed  approximately  2 miles  from 
threshold.  Readings  are  telemetered  or  relayed 
via  telephone  lines  to  central  facility,  along  with 
airport  centerfield  anemometer  reading.  A pilot  is 
notified  by  the  local  controller  if  the  vector  dif- 
ference between  the  distant  sensor  and  the  center- 
field  sensor  exceeds  some  threshold,  such  as  15  knots. 

Advantages:  Low/moderate  cost  and  maintenance;  uses  familiar 
and  relatively  simple  equipment.  Detects  outflow 
associated  with  gust  fronts  or  downburst  cells. 

Detects  passage  of  fronts  having  sharply  defined 
wind  shift  boundaries. 

Limitations : Does  not  provide  advance  warning  of  shear  events 
moving  into  the  corridor;  does  not  provide 
location  of  aircraft  encounter  with  frontal 
boundary  (slope  information  is  not  obtained); 
does  not  detect  synoptic  shear  above  airport 
environs . 
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Problems  Establish  whether  tradeoff  between  false 

Outstanding: alarms  caused  by  gusts  and  time  lag  associated 

with  smoothing  is  feasible;  establish  a threshold 
which  results  in  an  acceptable  miss/false  alarm 
tradeoff;  establish  whether  density  of  anemometers 
is  sufficient  to  prevent  misses  of  small  wind  shear 
event  s . 

c . Array  of  Anemometers  Around  Airpo rt 

Technique:  Anemometers  are  placed  as  in  a and  b above,  and 

additional  anemometers  are  placed  farther  out  and 
in  other  sectors.  Readings  are  relayed  by  telephone 
line  or  microwave  link  to  central  facility.  A pilot 
is  notified  if  wind  patterns  indicate  a shear  event 
moving  into  an  active  approach/departure  corridor. 
Pilots  are  alerted  as  to  ground  wind  speed  and  direc- 
tion at  runway  threshold.  Departing  aircraft  would 
be  advised  if  anemometer  readings  at  opposing  ends 
of  the  runway  indicated  significantly  different  winds, 
or  if  a shear  event  were  located  at  or  moving  into  the 
departure  corridor. 

Advantages:  'loderate  cost  and  maintenance;  uses  familiar  and 
relatively  simple  sensor  equipment.  Detects 
outflow  associated  with  gust  fronts  or  downburst 
cells.  Detects  passage  of  frontal  boundaries. 

Provides  advance  warning. 

Limitations : Complex  network  of  sensors  may  prove  costlv  to 

maintain;  real  estate  and  availability  of  vandalproof 
structures  may  prove  difficult  to  acquire.  Does  not 
provide  location  of  aircraft  encounter  with  frontal 
boundary;  does  not  detect  synoptic  shears  above 
airport  environs. 
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Problems  Establish  whether  resulting  wind  patterns  provide 
Outstanding:  signatures  that  can  be  identified  bv  a central 
processor;  establish  a threshold  and  averaeing 
time  which  results  in  an  acceptable  miss/false 
alarm  tradeoff. 

d . Array  of  Pressure  Jump  Sensors 

Technique:  Pressure  jump  sensors  are  placed  in  a grid  surrounding 
the  airport.  Readings  are  relayed  by  telephone  line 
to  central  facility.  Pressure  jumps  indicate  onset 
of  thunderstorm  gust  front,  and  time  sequences  of  jumps 
are  analyzed  to  estimate  time  of  arrival  of  gust  front 
at  operating  corridors.  A pilot  is  alerted  if  a 
shear  event  is  located  in,  or  is  moving  into  active 
corridor. 

Advantages:  Advance  warning  provided  for  thunderstorm  gust 

fronts.  Low/moderate  cost,  using  relatively  simple 
sensors;  can  be  mounted  on  telephone  poles. 

Limitations : Does  not  provide  an  all-clear  indication  when  shear 
event  has  passed.  Does  not  measure  winds  directly. 

Does  not  detect  frontal  passage  or  any  synoptic  shear 
events . 

Problems  Establish  that  pressure  jumps  are  associated  wi'th 
Outstanding:  most  hazardous  wind  shears;  establish  that  false 
alarms  (e.g.,  caused  by  gravity  waves)  are  suf- 
ficiently infrequent  or  are  recognizable;  establish 
means  of  estimating  shear  intensity. 

e . Acoustic  Sodar 

Technique:  The  acoustic  sodar  is  installed  beside  a shelter  on 

the  airport  surface.  The  sodar  points  upward,  and 
measures  the  two  horizontal  wind  components  at  heights 
from  100-1500  feet  (30-450m).  Vertical  wind  shears, 
winds  aloft,  and  surface  winds  are  measured  and  made 
available  to  controllers  and  pilots. 

I 
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Advantages:  Detects  synoptic  shears  and  locates  altitudes  of 
aircraft  encounters  with  vertical  shear  events. 

Detects  passage  of  fronts  having  sharply  defined  wind 
shift  zones.  Moderately  high  costs. 

Limitations: Susceptible  to  noise  from  aircraft,  rain  patter. 

Requires  expensive  bunker  installations  to  obtain 
data  above  500  feet  (150  m) . Does  not  detect  gust 
fronts  or  downburst  cells;  does  not  detect  cold  fronts 
before  they  arrive  at  airport. 

Problems  Establish  operational  feasibility  in  noisy  airport 

Outstanding:  environment,  in  rain,  and  in  conditions  which  exist 
when  nocturnal  boundary  layer  shears  are  present. 
Assess  reliability  and  maintenance  costs. 

f . FM/CW  Microwave  Radar 

Technique:  The  radar  is  installed  beside  a shelter  on  the  airport 

surface.  The  radar  points  upward,*  and  measures  the 
three  wind  components  at  heights  from  50-10,000  feet 
(15  m - 1 km).  Vertical  wind  shears,  winds  aloft, 
and  surface  winds  are  measured  and  made  available 
to  controllers  and  pilots. 

Advantages:  Detects  synoptic  shears  and  locates  altitudes  of 
aircraft  encounters  with  vertical  shear  events. 

Detects  passage  of  fronts  having  sharply  defined 
wind  shift  zones. 

Limitations : Moderately  high  cost,  complexity,  antennas  are 

bulky,  require  some  ground  installation.  Does  not 
detect  gust  fronts  or  downburst  cells;  does  not 
detect  cold  fronts  before  they  arrive  at  airport. 


Actually  the  radar  is  pointed  about  30°  from  the  zenith.  Clutter 
prevents  the  antenna  from  pointing  up  the  glide  slope;  thus  its 
use  is  limited  to  synoptic  measurements. 
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Problems  Establish  whether  returns  are  obtained  in  noctur- 
Outstanding; nal  boundary  layer  atmospheric  conditions,  i.e., 
establish  all-weather  potential.  Assess  relia- 
bility and  maintenance  costs. 


g . CW  Laser  Radar 

Technique:  The  laser  is  placed  in  a shelter  on  the  airport 

surface.  The  laser  beam  points  upward,  and  measures 
the  three  wind  components  at  heights  from  50-1000 
feet  (15  m - 300  m)  above  the  surface.  Vertical  wind 
shears,  winds  aloft,  and  surface  winds  are  measured 
and  made  available  to  controllers  and  pilots. 

Advantages:  All-weather  operation.  Detects  all  synoptic  shears 
and  locates  altitudes  of  Aircraft  encounter 
with  vertical  shear  events.  Detects  passage  of 
fronts  having  sharply  defined  wind  shift  zones. 
Compact  scanning  and  sensor  equipment. 

Limitations : Does  not  detect  gust  fronts  or  downburst  cells; 

does  not  detect  cold  fronts  before  they  arrive  at 
airport.  Moderately  high  cost,  complexity. 

Problems  Establish  turn-key  operational  capability,  elimi- 

Outstanding: nating  need  for  frequent  cryogenic  and  laser  gas 

refills.  Assess  reliability  and  maintenance  costs. 

h.  Pulsed  Doppler  Microwave  Radar 

Technique:  A highly  sensitive  Doppler  radar  is  installed  on 

the  airport  surface,  and  scanned  near  the  horizon. 
Returns  from  precipitation  or  index-of-refraction 
variations  in  the  atmosphere  permit  a line-of-sight 
wind  profile  to  be  obtained.  Horizontal  shears  along 
the  radials  are  detected,  and  wind  shear  events  can 
be  tracked.  Controllers  and  pilots  are  alerted  if 
windshear  events  are  located  in  or  are  approaching 
active  runway  corridors. 
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Advantages:  System  range  in  precipitation  is  more  than  adecuiatc; 

update  rate  is  likewise  more  than  adequate.  The 
shear  component  measured  corresponds  essentially  to 
headwind  gradients  encountered  by  the  aircraft. 
Provides  advance  warning  for  thunderstorm  gust  fronts 
and  downburst  cells.  Air  traffic  control  equipment 
and  maintenance  personnel  shelters  can  be  shared. 

L im i ta t i ons : U t i 1 i ze s large  (10-15  feet,  or  3-5  meters)  antennas 
requiring  tower  installation;  moderately  high  cost, 
moderate  complexity.  Vertical  beam  size  (somewhat 
greater  than  one  degree)  results  in  a limitation  on 
pointing  near  horizon  (caused  by  clutter)  and  results 
in  masking  of  shears  associated  with  winds  having 
significant  vertical  gradients.  Does  not  directly 
measure  downdrafts. 

Problems  Establish  capability  of  obtaining  clear  air 

Outstanding:  returns  in  thunderstorm  gust  front  and  downburst 
cells  where  precipitation  is  nonexistent. 

Establish  reasonable  cost/size/  performance 
tradeoffs  to  assess  ultimate  costs.  Establish 
signature  recognition  capability  and  false 
alarm  (e.g.,  due  to  gusts)  rejection  capabil- 
i ty . 

i . Pulsed  Doppler  Laser  Radar 

Technique:  A narrow-beam  laser  radar  is  placed  in  a shelter  on 

the  airport  surface  or  on  top  of  a building,  and 
scanned  close  to  the  horizon.  Returns  from  particu- 
late aerosols,  fog,  and  rain  permit  a 1 ine-of - s i ght 
wind  profile  to  be  obtained.  Horizontal  shears  along 
the  radials  are  detected,  and  wind  shear  events  are 
tracked.  Controllers  and  pilots  are  alerted  if  wind 
shear  events  are  located  in  or  are  approaching  active 
runway  corridors. 


3-36 


Advantages:  Provides  advance  warning  for  thunderstorm  gust  fronts 
and  downburst  cells.  The  shear  component  measured 
corresponds  essentially  to  headwind  gradients  en- 
countered by  aircraft.  Compact  equipment  can  be 
placed  in  a variety  of  locations  without  installation 
costs  associated  with  heavy  equipment.  Beam  can  be 
raised  in  elevation  to  provide  synoptic  shear  measure- 
ments. Beam  is  extremely  narrow,  and  can  be  pointed 
very  close  to  horizon  without  clutter  problems  or 
masking  of  non-horizontal  shears. 

Limi tations : Range  is  limited  in  heavy  rain,  fog,  and  snow. 

Does  not  directly  measure  downdrafts.  Does  not 
measure  winds  at  or  above  the  sensor.  Moderately 
high  cost,  complexity. 

Problems  I-stablish  signature  recognition  capability  and 
Outstanding: false  alarm  (e.g.,  due  to  gusts)  rejection 
capability;  establish  system  range  capabi- 
lity in  thunderstorm  environment.  Establish  turn- 
key operational  capability;  eliminate  need  for  fre- 
quent cryogenic  and  laser  gas  refills.  Assess 
reliability  and  maintenance  costs. 

Clearly,  many  of  the  limitations  listed  above  are  inherent: 
i.e.,  the  sensors  are  proposed  to  solve  only  one  problem.  Thus 
used  in  combinations,  they  cover  the  gamut  of  wind  shear  events. 
For  example,  at  Dulles  Airport  (Langweil,  1976)  a combination  of 
acoustic  sodar,  a small  pulsed  doppler  radar,  and  an  array  of 
pressure  jump  sensors  are  utilized  to  provide  advance  warning  of 
gust  fronts,  and  to  measure  synoptic  conditions. 

The  combination  of  a CW  laser  and  a pulsed  laser  complement 
each  other.  The  pulsed  laser  does  not  measure  winds  above  the 
airport;  although  the  elevation  angle  can  be  raised  to,  say,  6 
degrees,  and  synoptic  shears  can  be  detected  under  most  conditions 
under  conditions  of  heavy  fog  or  rain  the  penetration  would  be 


inadequate.  Raising  the  elevation  angle  higher  results  in  a 
loss  of  altitude  resolution.  Thus  a vertical-looking  CW  laser 
radar,  which  operates  in  all  weather  conditions,  would  measure 
the  vertical  wind  shears  in  fog  conditions  which  frequently 
accompany  warm  fronts,  where  a pulsed  laser  could  not.  The 
fortuitous  feature  of  the  laser  approach  is  that  most  of  the 
equipment  which  comprises  the  system  are  common:  scanner,  tele- 
scope, optics,  detector,  signal  processor,  and  data  processor. 
The  CK  laser  can  a] so  function  as  a master  oscillator  for  the 
pulsed  system.  The  result  is  that  the  capital  equipment  needed 
to  achieve  pulsed  operation  over  and  above  that  required  for  CW 
operation  costs  is  on  the  order  of  20  percent  or  less  of  the  CW 
system  cost.  Furthermore,  it  is  possible  to  switch  quickly 
between  CW  and  pulsed  modes  of  operation,  thus  achieving  a 
hybrid  CK/pulsed  operation  with  no  loss  of  performance  except 
the  time  lost  by  the  time-sharing  scheme. 

The  sections  which  follow  are  addressed  to  the  particular 
features,  problems,  and  characteristics  of  laser  radar  systems. 
The  system  requirements  developed  up  to  this  point  apply  to  any 
ground  sensor,  but  the  remaining  sections  develop  the  particular 
laser  design  constraints  these  requirements  entail. 


SENSORS  AND  PROPAGATION  CHARACTERISTICS 


4.1  BASIC  FEATURES  OF  THE  MEASUREMENT 

The  fundamental  property  of  the  atmosphere  that  is  exploited 
by  a laser  radar,  or  lidar,  is  the  presence  of  particles  which  are 
large  in  comparison  with  the  air  molecules,  and  move  essentially 
with  the  air  mass.  The  particles,  which  are  plentiful  except  in 
rare  dry,  clean  conditions,  are  dust  particles,  moisture  droplets, 
fog  droplets,  raindrops,  or  particulate  pollution,  frequently  all 
grouped  under  the  name  "aerosol."  The  laser  transmitter  directs 
a beam  toward  a narrow  volume  of  air;  some  radiation  is  scattered 
and  some  absorbed  by  each  aerosol  in  the  beam  path,  and  thus 
never  reaches  the  backscattering  volume;  this  causes  an  attenua- 
tion of  the  beam.  Some  of  the  radiation  reaching  the  backscat- 
tering volume  is  scattered  back  toward  the  laser  (backscatter) 
by  aerosols  in  the  volume,  and  is  again  attenuated  by  aerosols  in 
the  return  path  (see  Figure  4-1).  The  signal  received  is  detected 
by  heterodyning  with  a local  oscillator  whose  frequency  is  tightly 
controlled.  The  Doppler  shift  principle  dictates  that  the  fre- 
quency of  the  backscattered  signal  differs  from  that  of  the  trans- 
mitted signal  by  an  amount  proportional  to  the  velocity  component, 
parallel  to  the  beam,  of  the  particles  in  the  volume.  Thus  the 
measurement  of  this  frequency  difference  is  a measure  of  the 
velocity  of  the  air  mass  of  the  volume  in  the  direction  of  the  beam. 

It  is  informative  to  consider  the  effect  of  aerosol  densities 
on  the  received  signal  (see  Figure  4-2).  Suppose  the  volume  of 
interest  is  fixed  at  some  distance  from  the  laser,  and  the  density 
of  particles,  of  a given  size  distribution,  increases  from  zero: 
at  very  low  density,  attenuation  is  negligible,  but  backscatter 
is  very  small,  resulting  in  a small  signal;  as  the  density  in- 
creases, the  attenuation  remains  small,  backscatter  is  high, 
resulting  in  a large  signal;  eventually  a point  is  reached  where 
attenuation  reduces  the  signal  at  the  same  rate  that  backscatter 
increases  it,  and  the  signal  is  maximized;  beyond  that,  attenua- 
tion reduces  the  signal  much  faster  than  backscatter  increases 
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it,  resulting  in  a rapidly  decreasing  signal.  Some  fogs  and 
heavy  rain  can  cause  loss  of  signal  at  2-3  nautical  miles  range 
(4-6  kilometers)  with  present  day  lasers. 

There  are  several  sources  of  error  which  can  be  troublesome 
if  caution  is  not  exercised  in  data  handling:  birds  and  insects 
flying  through  the  volume  of  interest  can  cause  sudden  large 
signals  with  erroneous  velocities;  falling  rain  can  have  line- 
of-sight  components  that  differ  from  the  winds;  also,  turbulence 
can  cause  spatial  and  temporal  variations  in  the  wind. 

Carbon  dioxide  (CO^)  gas  lasers  operating  at  10.6  micron 
wavelengths  are  the  only  type  that  will  be  considered  in  this 
study.  This  is  chiefly  because  of  the  technology  available  only 
at  this  wavelength:  high  power,  stable  single-mode  operation, 
high  Uv’^ector  efficiency  and  convenient  video  bandwidths.  There 
are  propagat ional  advantages  as  well:  penetration  in  light  fog 
is  markedly  superior,  and  susceptibility  to  turbulence  is  signi- 
ficantly lower  than  at  smaller  wavelengths, 

4.2  PROPAGATION  PROPERTIES  OF  THE  ATMOSPHERE 

In  order  to  be  a viable  operational  sensor,  the  laser  system 
must  be  able  to  operate  under  a wide  range  of  conditions:  varying 
humidities,  aerosol  content,  haze,  fog,  and  precipitation.  These 
vary  geographically  and  seasonally.  Conditions  range  from  very 
dry,  clean  air  where  backscatter ing  is  small,  to  snowstorms,  where 
attenuation  is  extreme.  Fortunately,  aerosols  are  concentrated 
near  the  ground,  and  most  operational  environments  have  a signi- 
ficant aerosol  content  when  hazardous  wind  shears  are  likely  to 
be  present.  At  the  other  extreme  where  rain  limits  laser  range, 
shear-producing  thunderstorm  cells  have  highly  inhomogeneous  rain 
patterns;  the  gust  fronts  associated  with  the  cells  precede  the 
regions  of  heavy  precipitation  and  can  thus  be  detected  prior  to 
the  arrival  of  the  rain  cell. 
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There  are  several  physical  mechanisms  which  affect  the 
amplitude  of  the  received  signal:  (1)  attenuation,  by  aerosol 
scattering  and  absorption;  (2)  backscatter  by  the  aerosol  assem- 
blage; (3)  fluctuations  of  the  backscattered  signal;  (4)  fluctua- 
tions of  the  attenuation  by  scattering;  (5)  scintillation,  i.e., 
amplitude  fluctuations  caused  by  index-of-refraction  and  aerosol 
inhomogeneities  in  the  beam  path;  and  (7)  beam  spread,  caused  by 
scattering. 

Of  the  7 factors,  only  (l)-(3)  and  (6)  are  believed  to  be 
significant  in  the  present  application.  Fluctuations  of  the 
attenuation  are  small.  Scintillation  is  likewise  small,  even 
with  considerable  turbulence.  Angle-of -arrival  variations  are 
not  significant,  because  reciprocity  requires  the  transmitted 
and  received  paths  to  be  the  same.  Phase  front  distortion, 
however,  can  be  significant  in  heavy  turbulence,  since  it  reduces 
the  effective  receiving  aperture.  Beam  spread  by  scattering 
occurs  under  high  attenuations,  where  range  is  primarily  limited 
by  the  attenuation;  the  beam  spread  signal  loss  is  small  compared 
with  the  attenuation  loss. 

4.2.1  Attenuation 

F.ach  particle  between  the  transmitter  and  the  volume  of 
interest  scatters  some  energy  into  a 4tt  steradians  of  solid 
angle;  in  addition  some  energy  is  absorbed  by  the  particle  and 
converted  to  heat.  Particles  smaller  than  a wavelength  scatter 
relatively  weakly  (Rayleigh  scattering)  while  particles  larger 
than  a wavelength  scatter  an  amount  approximately  proportional 
to  the  cross-sectional  area  (Mie  scattering).  Thus  for  a given 
beam  volume  the  attenuation  depends  on  the  particle  size  distri- 
bution, the  number  of  particles,  the  wavelength,  and  the  particle 
indices  of  refraction. 

In  clear  air  molecular  absorption,  rather  than  aerosol 
absorption  or  scattering,  appears  to  be  the  limiting  factor. 

The  values  for  attenuation  published  by  McClatchey,  et  al . (1972) 

1 

indicate  that  the  attenuation  by  molecular  absorption  ranges  from 

1 
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0.63  dB/nm  (0.34  dB/km)  in  winter  to  2.9  dB/nm  (1.6  dB/km)  in 
.summer.  Clear  air  aerosols  would  only  give  about  0.07  dB/nm 
(0.04  dB.kmJ,  while  haze  by  itself  would  cause  about  0.4  dB/nm 
(0.2  dB/km). 

In  light  developing  fogs  and  mist,  the  attenuation  is  small, 
typically  2-3  dB/nm  (1-1.6  dB/km);  under  these  conditions  the 
wavelength  of  10.6  microns  has  significantly  greater  penetrating 
capability  than  optical  wavelengths.  The  primary  reason  for  this 
is  that  the  particles  are  primarily  of  a size  which  results  in 
weak  (Rayleigh)  scattering  at  infrared  and  strong  (Mie)  scatter- 
ing at  optical  frequencies.  In  denser  fogs,  the  average  drop 
size  also  increases,  and  attenuation  can  reach  extreme  values 
of  up  to  75  dB/nm  (40  dB/km).  For  a given  water  density  the 
attenuation  at  10.6  micron  wavelength  is  relatively  insensitive 
to  drop  size  distribution  (but  highly  sensitive  at  optical  wave- 
lengths); it  can  be  approximated  by: 

a = 0.5  p dB/km 

where  p is  the  water  density  ^.n  milligrams  per  cubic  meter. 

If  p is  given  in  pounds  per  cubic  foot,  a is  approximately; 

a = 6x10  ^p  dB/nm  . 

A typical  light  fog  has  a water  density  of  5 mg/m^  (3x10’^  lb/ 
Ct^)  while  a typical  dense  fog  may  reach  100  mg/m^  (6x10  ^ lb/ 
ft^). 

In  rain,  only  Mie  scattering  occurs;  the  attenuation  varies 
inversely  witli  the  drop  diameter  for  a given  water  density.  Drop 
sizes  typically  range  between  0.2  mm  (0.008  in)  and  6 mm  (0.2 
in)  in  diameter.  Rainfalls  observed  in  nature  have  a wide 
spread  of  drop  sizes  for  any  rainfall  rate,  but  they  can  be 
characterized  by  the  following  behavior:  the  most  prevalent 
drop  size  increases  with  rain  rate.  This  is  shown  in  Figure  4-3, 
which  is  taken  from  Kerr  (1951),  based  on  measurements  by  Laws 
and  Parsons  (1943).  Using  this  model,  Chu  and  Hogg  (1968) 
calculated  the  attenuation  as  a function  of  rainfall  rate;  this 
is  shown  in  Figure  4-4.  It  is  evident  that  heavy  rainfall  rates 
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jircatly  reduce  the  signal  strength. 


4.2.2  Backscatter 

liach  particle  in  the  volume  of  interest  scatters  some  energy 
back  toward  the  telescope.  As  in  the  case  of  attenuation,  weak 
(Rayleigh)  scattering  occurs  in  clear  air  and  haze,  due  to  the 
existence  of  a few  aerosols  which  are  usually  present.  Fog  and 
rain  result  in  an  enhanced  backscattered  signal,  but  since 
attenuation  also  increases,  the  net  result  is  a reduced  signal 
at  long  range. 

The  backscatter ing  coefficient  6(11)  varies  by  several  orders 

of  magnitude,  depending  on  conditions.  In  clear  air  and  haze 

- 9 - 6 2 3* 

B(ti)  generally  varies  from  10  to  10  m /sr/m  . In  fog,  6(Tt) 

"6  ”52  3 

varies  from  about  10  to  4x10  m'^/sr/m'^;  in  rain,  from  about 

j 5 2 3 

10  to  2x10  m^/sr/m^.  The  backscattering  coefficient  expresses 

the  equivalent  cross-sectional  area,  presented  by  a cubic  meter 
of  scattering  particles  to  the  incident  beam,  that  would  scatter 
the  same  fraction  of  energy  per  unit  of  solid  angle  back  toward 
the  transmitter. 

There  is  obviously  a connection  between  the  backscatter 
coefficient  and  attenuation  coefficient  (a),  since  they  exhibit 
the  same  dependence  on  dielectric  constant  and  drop  diameter. 

However,  the  scattering  diagram  of  a drop  varies  considerably  as 
a function  of  drop  size.  But,  since  these  variations  are  due  to 
resonances,  the  fact  that  drops  of  varying  sizes  are  present  in 
a. volume  has  the  effect  of  averaging  the  resonances,  and  a 
reasonable  estimate  at  a relation  between  the  two  averaged 
coefficients  can  be  made  for  water  droplets  larger  than  a wave 
length.  This  is  done  in  Appendix  B.  The  result  is: 
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where  a is  the  attenuation  in  nepers  per  meter.  (a  in  nepers 
per  meter  is  equal  to  a in  dB/km  divided  by  4343,  and  to  a in 
dB/nm  divided  by  8035.) 

The  relation  between  a and  B(tt)  used  here  holds  for  lari;e 
droplets,  but  the  assumptions  break  down  for  droplet  diameters 
of  loss  than  5-10  microns.  Nonetheless,  it  is  still  a good 
working  figure  down  to  about  2 microns;  below  that,  attenuation 
is  due  mostly  to  absorption,  not  scattering.  Backscattering  is 
strongly  dependent  on  drop  diameter,  while  absorption  is  not, 
so  that  the  ratio  of  BC^r)  to  a decreases  for  smaller  droplets. 
This  means  that  for  attenuations  less  than  2 dB/nm  (1  dB/km), 
backscattering  estimates  are  too  high,  the  actual  values  being 
highly  dependent  on  the  drop  sizes  present:  the  smaller  the 
droplets,  the  smaller  the  backscattering  coefficient. 

At  the  one  extreme  of  heavy  rain  or  dense  fog,  the  signal 
at  any  significant  range  is  limited  by  attenuation.  At  the  other 
extreme,  the  signal  is  limited  by  the  lack  of  backscatter.  The 
latter  situation  does  not  appear  to  be  a serious  problem,  because 
most  airports  are  near  urban  centers,  where  particulate  pollution 
is  present  or  near  salt  water  where  salt  particles  act  as  hvdro- 
scopic  nuclei.  In  fact,  even  in  relatively  clean  air,  there  are 
generally  a sufficient  number  of  aerosols  between  1-20  microns  in 
diameter  to  provide  a detectable  signal.  An  extreme  example  was 
observed  in  a recent  measurement  program  (Brashears,  et  al , 1977) 
utilizing  a CW  laser.  After  a cold  front  passed  through  on  a 
dry  winter  day  at  Table  Mountain  in  Boulder,  Colorado,  a notice- 
able drop  in  signal  occurred:  the  backscatter  coefficient  was 
estimated  to  be  about  10  m^.  Even  here,  the  signal  was 

sufficient  to  establish  the  wind  velocity.  These  relatively  rare 
conditions  are  not  believed  to  be  associated  with  hazardous  wind 
shears . 
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A phcnonicMion  w'v  ich  is  encountered  by  a focused  laser  is 
that  for  rain  the  backscatter ing  volume  may  be  so  small  that  on 
the  average  less  than  one  raindrop  is  present  at  a time.  In 
this  case  the  backset ttering  is  provided  by  aerosols;  however, 
there  are  still  many  particles  in  the  path  to  cause  attenuation. 
This  phenomenon  does  not  occur  in  fog,  where  the  particle  density 
is  much  higher;  nor  does  it  occur  for  pulsed  lasers  (focused  at 
infinity)  where  the  volume  is  orders  of  magnitude  larger. 

4.2.3  Fluctuations  of  the  Backseat  ered  Signal 

The  backscattered  signal  is  the  sum  of  the  returns  from 
particles  in  the  volume  of  interest.  On  the  scale  of  the  inter- 
particle  spacing,  the  locations  of  the  various  particles  are 
random  and  independent.  Thus  the  phases  of  their  returns  can  be 
considered  to  be  random  and  independent.  This  gives  rise  to  a 
combined  signal  exhibiting  Rayleigh  statistics.  If  all  the 
particles  in  the  volume  are  moving  with  the  same  velocity,  the 
frequency  of  the  combined  signal  will  be  Doppler-shifted  from 
the  transmitted  frequency  by  a fixed  amount,  and  the  amplitude 
will  be  random,  with  a Rayleigh  probability  distribution.  If 
w is  the  rms  value  of  the  intensity  of  the  return  from  each 
particle,  the  mean  intensity  is  given  by  ttNw/4,  where  N is  the 
number  of  particles.  The  standard  deviation  of  the  intensity  is 
given  by  Nw . Thus,  the  standard  deviation  approximates  the  mean. 
This  means  that  the  signal  fluctuates  quite  widely.  The  probabil- 
ity that  the  received  signal  (s)  exceeds  some  intensity  value 
(.SJ  is  given  by 

P (s  > S)  = 1 - exp  (-S^tt/4M) 

where  M is  the  mean  intensity;  this  is  shown  in  Figure  4-5. 

Irom  this  it  can  be  seen  that  5 percent  of  the  returns  will  either 
exceed  the  mean  by  6 dB  or  be  down  12  dB  or  more  below  the  mean. 

The  question  of  coherent  detection  naturally  arises:  is 
the  coherence  time  of  the  phenomenon  long  enough  that  waveforms 
from  successive  pulses  could  be  combined  to  reduce  the  effect  of 
noise?  The  answer  appears  to  be  negative,  since  in  the  presence 
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i 

1 

1 

of  a shear,  different  portions  of  the  volume  have  a relative 
motion  which  destroys  coherence  on  a time  scale  of  microseconds. 

The  shear  also  causes  a spectral  spread  (another  manifestation  of 
the  same  thing) . 

This  does  not  say  that  averaging  is  not  helpful;  incoherent 
averaging  is  useful  because  while  the  noise  mean  does  not  improve, 
the  noise  variance  is  reduced.  This  allows  postintegration  j 

thresholding  to  be  performed  prior  to  data  processing,  where 

parameter  estimation  processing  takes  place.  Also,  the  processor  > 

can  be  designed  to  ignore  samples  where  the  signal-to-noise  ratio  1 

(SNR)  drops  below  some  threshold  (e.g.,  3dB) , so  that  signal  ; 

dropouts  do  not  degrade  the  averaged  SNR.  | . 

4.2.4  Fluctuations  Due  to  Beam  Path  Disturbances 

As  a laser  beam  travels  through  the  atmosphere,  there  are 
a number  of  mechanisms  which  alter  its  shape,  direction  and 
intensity.  Depending  on  the  cause  of  these  disturbances, 
different  characteristics  are  affected  at  different  wavelengths. 

At  wavelengths  near  10.6  microns,  turbulence  effects,  which 
cause  i ndex-of- ref raction  variations,  are  much  smaller  than  at 
optical  wavelengths;  they  must  be  considered,  however.  The 
attenuation  and  changes  in  attenuation  caused  by  turbulence  are 
not  significant.  Phase  front  disturbances,  on  the  other  hand, 
can  limit  performance.  Small  eddies  encountered  by  the  propagat 
ing  beam  cause  the  beam  to  spread,  and  cause  phase  front  pertur- 
bations, while  large  eddies  actually  deflect  the  beam  (Fante, 

1075).  For  a heterodyned  laser  beam  the  beam  deflection  (also 
called  beam  wander  or  spot  dance)  is  not  a problem,  since  the 
effect  is  merely  to  change  the  volume  of  scatterers  being  illumi 
nated:  by  reciprocity,  the  return  signal  traverses  the  same  path 
back  to  the  telescope.  Beam  spread  caused  by  turbulence  is  also 
small  (see  equation  (37)  of  Fante,  1975).  The  phase  front  per- 
turbations have  the  effect  of  limiting  the  dimension  across  the 
phase  front  over  which  coherence  is  maintained.  Stated  another 
way,  for  a given  rms  phase  ripple  increasing  telescope  diameter 
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will  not  be  effective  (floldstein.  Miles,  and  Chabot,  1965  and 
Strohbchn,  1968). 


Thus  the  effect  of  turbulence  can  be  summed  up  as  limiting 
the  effective  diameter  of  the  collecting  optics.  This  elTective 
diameter  is  given  by  (Goldstein  et  al . , 1965): 

5/5 


n / 0.0588  ) 


meters 


where  X is  the  wavelength  in  meters,  is  a measure  of  the 
level  of  turbulence  having  the  units  of  -2/3  meters,  and  K is 
the  path  length  (two-way)  in  meters. 


For  typical  values  of  turbulence  near  the  ground,  is 

about  4x10  m which  yields  an  effective  diameter  of  0.76  m 

(2.5  ft)  for  26  km  (14  nm)  of  path  length.  For  strong  turbulence, 
is  about  10  m which  yields  an  effective  diameter  of 

0.11m  (0.36  ft).  In  order  that  1 foot  (0.3m)  optics  be  effec- 
tive at  ranges  of  7 nm  (26  km  two-way  path),  the  average  value 
of  would  have  to  be  less  than  2x10  Whether  this 

poses  a serious  limitation  is  not  known. 


Turbulence  which  causes  the  index-of-refraction  variations 
is  associated  with  wind  shears  and  thermal  convective  heating 
from  the  ground  (Strohbehn,  1968).  The  latter  causes  turbulence 
which  is  strongest  in  the  lowest  100  feet  (30  m) , and  thus  is 
less  of  a problem  for  a laser  beam  which  will  be  looking  up  the 

glide  slope.  However,  whether  values  of  in  excess  of 

- 1 5 - 2/3  " 

2x10  m are  encountered  in  thunderstorm  outflows  is  not 

known . 

Particle  scattering  primarily  causes  attenuation.  Systematic 
density  variations  in  the  transverse  direction  could  in  principle 
cause  beam  bending,  but  there  is  no  reason  to  believe  that 
particle  densities  exhibit  such  variations.  Forward  scattering 
could  conceivably  cause  beam  spread,  and  this  has  been  observed 
in  measurements  in  rain  and  snow  (Chu  and  Hogg,  1968).  However, 
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since  the  attenuations  calculated  for  particle  scattering  already 
account  for  the  energy  not  available  to  the  barkscattering  volume, 
it  does  not  appear  that  the  beam  spread  represents  an  additional 
loss  over  and  above  the  attenuation  value. 

Phase  fluctuations  caused  by  particle  scattering  are  believed 
to  be  small,  but  there  are  no  rigorous  calculations  or  experiments 
which  ver i fy  this . 

4.3  CW  LASERS  IN  THE  VAD  SCAN  MODE 
4.3.1  Basic  Measurement  Technique 

A CW  heterodyned  laser  can  obtain  wind  profiles  in  a volume 
above  the  laser  by  employing  a VAD  (Velocity  Azimuth  Display) 
scan,  similar  to  that  described  by  Lhermitte  and  Atlas  (1961) 
lor  radar.  In  this  mode  the  range  is  varied  by  focusing  the  beam 
at  the  range  of  interest,  and  scanning  the  beam  in  azimuth,  while 
holding,  the  angle  from  the  zenith  constant  (typically  30  degrees) 
(see  Figure  4-6).  If  the  horizontal  wind  remains  constant  during 
the  time  of  a complete  rotation,  the  line-of-  sight  Doppler 
velocity  component  recorded  will  be  a sine  wave  whose  amplitude 
is  a measure  of  horizontal  wind  direction,  and  whose  average 
value  is  a measure  of  vertical  draft  speed.  By  refocusing  the 
telescope  to  a series  of  heights,  a profile  of  the  three  wind 
components  can  be  mapped  out. 

Deviations  from  a sine  wave  denote  nonuniformity  of  the 
wind  field,  or  changes  in  the  wind  field  with  time.  Turbulence 
would  cause  both  of  these  variations,  so  that  the  nonuniformity 
can  potentially  be  used  as  a measure  of  turbulence.  This  is 
discussed  further  in  Section  4.3.6. 

Laser  systems  have  been  fabricated  and  tested  which  use 
this  principle  (Brashears  and  Hallock,  1976).  They  have  been 
successfully  employed  and  tested  against  tower  anemometer 
measurements.  As  a result,  the  requirements  of  the  following 
sections  arc  based  on  the  observed  behavior  of  an  operating 
sy  s t em . 
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FIGURE  4-6.  VELOCITY-AZIMUTH-DISPLAY  (VAD) 
SCAN  CONCEPT 
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4.3.2  .Scanning  Tcchniquc.s 


Ihere  arc  two  scanning  modes  which  may  be  employed:  (1) 

King  .Scan,  which  is  a rapid  azimuth  scan  (e.g.,  0.5  revolution 
per  second)  combined  with  a stepped  range  scan  (e.g.,  at  50 
foot  (15  mj  increments  ) (see  Figure  4-7);  or  (2)  Finger  Scan, 
which  is  a rapid  range  scan  (e.g.,  0.3  Hz)  combined  with  a slow 
azimuth  scan  (e.g.,  2 revolutions  per  minute)  (see  Figure  4-8), 
Ihere  is  a wide  latitude  of  choices,  because  for  synoptic 
phenomena,  changes  take  place  slowly,  and  there  is  no  need  for 
a rapid  update  rate.  A CW  laser  sensor  operating  full-time 
would  utilize  extensive  averaging,  and  take  much  longer  than  the 
40  second  update  rate  used  to  derive  the  parenthetical  nunibers 
listed  above. 

The  reason  for  minimizing  the  total  measurement  period  is 
in  anticipation  of  the  combined  CW/pulsed  operation;  there  the 
need  for  rapid  updating  of  information  for  the  pulsed  mode  of 
operation  is  such  that  as  little  time  as  necessary,  for  a 
reliable  synoptic  measurement,  should  be  allocated  to  the  CW 
mode  of  operation.  A measurement  of  40  seconds  duration  every 
10-15  minutes  is  probably  adequate. 

4.3.3  Altitude  Resolution 

The  rationale  for  an  altitude  resolution  has  already  been 
developed  in  Section  3.3,  where  a 50  foot  (15  m)  altitude 
resolution  is  recommended  to  properly  characterize  frontal 
shears  (see  especially  Figure  3-9). 

4.3.4  Accuracy 

In  section  3.4  the  recommended  system  accuracy  was 
determined  to  be  1-2  kphf  (0.017-0.034  m/sec/m).  Since  the 
measurement  technique  yields  the  magnitude  and  phase  of  the 
horizontal  winds  (plus  net  vertical  flow),  a polar- to-rectangular 
transformation  is  necessary,  where  the  rectangular  coordinate 
axes  are  aligned  with  the  runway  in  order  to  provide  headwind 
(and  crosswind)  profiles.  From  the  headwind  profile  successive 
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differences  between  altitudes  at  50  foot  (15  m)  intervals  would 
give  the  shear: 

WS(h)  * 2(HW(h-25)  - HW(h+25)^kphf 

where  h is  the  height  AGL,  height  units  are  in  feet,  and  head- 
winds are  in  knots;  a positive  shear  indicates  an  increasing 
headwind  for  a descending  aircraft,  and  vice  versa,  (An 
ascending  aircraft  would  experience  a decreasing  headwind  in 
the  same  situation.) 

In  order  to  achieve  2 kphf  of  accuracy  using  this  formula, 
the  winds  would  have  to  be  known  to  an  accuracy  of  0.71  knot 
(0.36  m/sec),  assuming  the  two  wind  values  to  be  independent. 
I’aradoxical  ly , if  wind  profile  values  were  only  available  at 
100  foot  (30  m)  intervals,  the  wind  shear  would  be  given  by: 

WS(h)  = HW(h-50)  - llW(h-*-50)  kphf 

and  2 kphf  accuracy  would  be  achieved  with  an  accuracy  of  1.4 
knots  (0.7  m/sec).  This  would  suggest,  at  first  glance,  tliat 
accuracy  could  be  improved  by  ignoring  alternate  measurements. 
The  key  to  this  paradox  lies  in  the  assumption  of  independence 
of  wind  values  at  adjacent  altitudes.  In  fact,  the  accuracy 
can  be  significantly  improved  by  proper  smoothing  of  the  two- 
dimensional  data.  This  is  di.scussed  in  section  4.3.7.  As  a 
guess,  the  accuracy  of  2 kphf  can  probably  be  achieved  with 
a wind  profile  measurement  accuracy  of  about  2 knots  (1  m/sec), 

4.3.5  Signal  Processing 

The  signal  processing  operates  on  the  raw  signal,  and 
produces  a frequency  domain  set  of  descriptors:  either  the 
entire  spectrum,  or  a set  of  parameters,  e.g,,  the  strongest 
frequency,  the  amplitude  at  the  strongest  frequency,  and  the 
highest  and  lowest  frequencies  of  the  spectral  bandwidth.  The 
processing  can  be  thought  of  as  a Fourier  transform  operation. 
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with  or  without  interpretation.  In  addition,  signal  process- 
ing may  include  rejection  of  low-amplitude  frequencies. 

The  sample  period  is  quite  arbitrary,  but  is  typically 
tens  of  milliseconds.  The  spectra  which  result  from  the 
processing  by  a spectrum  analyzer  typically  look  like  Figure 
4-9a  and  b:  the  first  spectrum  is  for  a steady  wind,  where  only 
one  wind  velocity  is  present;  the  second  spectrum  is  broadened, 
which  could  be  caused  by  a shear  over  the  volume,  or  by 
turbulence,  where  several  wind  velocities  are  present. 

For  one-foot  (0.3  m)  optics,  the  focal  volume  is  approx- 
imately a cylinder  whose  diameter  varies  from  0.04  inch  (1  mm) 
at  a 50  foot  (15  m)  range,  to  0.5  inch  (12  mm)  at  1000  feet 
(300  m) ; its  length  varies  simultaneously  from  0.6  foot  (0.2  m) 
to  200  feet  (60  m) , respectively.  Thus  near  the  ground  the 
focal  volume  is  so  small  that  individual  spectra  are  quite 
narrow,  and  turbulence  is  seen  rather  as  a sample- to- sample 
variation.  Higher  up,  however,  the  focal  volume  is  large 
enough  to  incorporate  a large  number  and  size  of  eddies,  and 
turbulence  will  be  seen  as  a spectral  broadening,  indicating 
that  there  are  wind  speed  variations  within  the  volume. 

The  characteristic  sine  wave  of  the  VAD  scan  is  a plot 
of  the  frequency  (proportional  to  the  wind  velocity)  versus  the 
scan  angle.  Thus  one  of  the  problems  is  to  choose  the  frequency 
parameter  that  best  characterizes  the  average  frequency  of  a 
given  spectrum.  Some  candidate  parameters  are  the  following: 

a.  Frequency  at  the  peak  amplitude  - This  is  simplest, 
but  noise  spikes  cause  significant  errors. 

b . Average  of  minimum  and  maximum  frequencies  of  the 
spectral  band  - This  is  better,  but  any  skew  tends  to 
be  accentuated  where  particulate  concentration  is  high. 

c . Average  of  minimum  and  maximum  frequencies  of  the 
spectral  band,  measured  at  some  fraction  of  the  peak 
amplitude  - This  is  a more  consistent  measurement,  but 
somewhat  complicated  to  implement. 
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FIGURE  4-9.  DOPPLER  WIND  SPECTRA 
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d.  Centroid  of  spectral  band  - Probably  the  best  measure- 
ment, but  complicated. 

It  should  also  be  pointed  out  that  the  difference  between 
signal  processing  and  data  processing  is  not  precise:  parameter 
determination  may  well  be  implemented  in  software. 

If  signal  processing  is  performed  in  hardware,  it  must  be 
able  to  handle  unusual  situations  like  insect  replies,  as  well 
as  returns  from  falling  rain.  Both  of  these  cause  "double- 
peaking"  in  the  spectrum,  since  there  are  two  or  more  distinct 
velocities  present  in  the  return.  Some  form  of  outlier  rejection 
scheme  in  the  data  processing  could  help  in  the  case  of  insects, 
while  choosing  the  slower  of  two  or  more  sustained  peaks  would 
result  in  the  correct  characterization  of  the  wind  in  the  presence 
of  rain. 

4.3.6  Data  Processing 

In  both  scan  modes  (Ring  Scan  and  Finger  Scan)  the  result 
is  a collection  of  Doppler  velocities  at  a number  of  azimuth 
angles  for  each  range,  or  height.  Figure  4-10  shows  a measure- 
ment example  (Brashcars,  et  al,  1976):  the  missing  values 
between  *2  knots  are  the  result  of  preprocessing  thresholding. 

It  is  evident  that  several  different  methods  could  be  used  for 
estimating  the  wind  speed  and  direction  from  the  data.  Some 
of  these  are  discussed  below: 

a.  Peak  Algorithm:  take  the  largest  value  of  wind  speed, 
and  use  that  value  as  the  wind  speed,  and  the  associated 
phase  as  wind  direction.  This  has  the  result  of  always 
estimating  high  in  wind  speed,  and  is  not  reliable  in 
direction. 

b.  Filtered  Peak  Algorithm:  use  some  averaging  or  filtering 
technique  to  reduce  the  effect  of  noise,  and  proceed  as 
in  "a"  above.  This  reduces  the  effect  of  noise,  but 
still  gives  errors  due  to  turbulence  and  gusts,  which 
cause  systematic  deviations  from  a sine  wave. 
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c.  Fourier  Transform  Technique:  recover  the  fundamental 
frequency  sine  wave  by  a digital  Fourier  transform  cal- 
culation, and  use  the  resulting  sine  wave  to  estimate 
wind  speed  and  direction.  This  takes  some  computation, 
but  provides  the  best  estimate. 

Algorithms  b and  c above  can  be  improved  by  some  form  of 
outlier  rejection  scheme,  whereby  data  which  vary  from  their 
neighbors  by  some  amount  are  rejected 

There  are  some  differences  between  the  two  scanning  modes 
that  favor  the  Finger  Scan  mode.  For  one  thing,  turbulent  eddies 
whose  characteristic  lifetimes  are  longer  than  the  Ring  Scan 
period  (e.g.,  2 seconds),  but  are  less  than  the  total  measurement 
period  (e.g.,  40  seconds)  will  tend  to  be  averaged  out  in  the 
linger  Scan  mode,  while  they  will  tend  to  bias  the  results  in 
the  Ring  Scan  mode.  For  another,  in  the  Finger  Scan  mode,  samples 
are  taken  at  more  altitudes  (and  less  azimuth  angles);  equivalent 
accuracy  and  resolution  with  the  Ring  Scan  mode  can  be  obtained 
by  averaging  the  samples  within  the  altitude  resolution  cell. 
Alternatively,  the  additional  altitude  samples  can  be  used  to 
obtain  a higher  altitude  resolution,  albeit  at  a reduced  accuracy. 
This  feature  may  be  useful  in  assessing  the  intensity  and  extent 
of  shear  zones. 

Once  a profile  is  obtained,  the  result  can  be  further 
scrutinized  to  sec  if  it  is  physically  realizable.  For  example, 
synoptic  profiles  approach  logarithmic  in  the  lowest  100  feet 
(30  m) ; a sequence  of  >10,  ♦lO,  -8,  >10  knots  (+5,  +5,  -4,  *5 
m/scej  at  50  foot  (15  m)  altitude  increments  would  not  be 
physically  realizable.  As  more  shear  data  are  obtained,  other 
physically  realizable  constraints  can  be  brought  in  to  improve 
confidence  in  predicted  shears  and  reduce  false  alarms  due  to 
spurious  measurements. 
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4.3.7  CW  Laser  Hazard  Limitations 


Laser  hazards  potentially  exist  within  the  shelter  and 
outside.  The  hazards  within  the  shelter  can  be  minimized  by 
proper  enclosures  and  precautions.  Outside  the  shelter,  there 
are  the  questions  of  hazards  to  maintenance  personnel  on  the 
roof  of  the  shelter  or  in  the  area  of  the  shelter,  people  in 
the  airport  buildings,  and  passengers  and  crew  in  aircraft. 

Considering  only  maximum  intensity,  there  are  regions 

where  the  beam  is  focused  (the  focal  volume)  where  intensities 

are  quite  high.  For  example,  the  worst  case  is  at  minimum 

range  (50  feet,  or  15  m) , where  the  beam  is  focused  down  to  a 

volume  1 mm  (0.04  in)  in  diameter.  There  a 10  watt  laser  would 

2 

result  in  a peak  intensity  of  1270  watts/cm  *,  which  is 
definitely  dangerous  to  skin  and  eyes;  however,  a maintenance 
man  on  the  roof  of  the  shelter  would  experience  at  most  only 
0.02  watt/cm^,  a safe. amount.  Furthermore,  if  the  scanner 
were  in  motion,  the  amount  of  time  that  a spot  would  be  illumi- 
nated would  be  small,  which  is  not  as  dangerous  as  continued 
irradiation. 

The  Food  and  Drug  Administration  has  issued  standards  on 
laser  radiation  (FDA,  1975),  which  define  three  classes  of 
at  CO^  wavelengths  (Class  li  does  not  apply  to  CO2  wavelengths): 

Class  I:  For  a scanned  laser,  the  radiant  energy  through 

a 7 mm  diameter  aperture  shall  not  exceed  4.4x10  ^ t 

joules  where  t is  the  duration  of  the  irradiation  in  seconds. 

For  a fixed  beam  laser,  the  radiant  power  through  an  80  mm 

- 4 

diameter  aperture  shall  not  exceed  7.9x10  watts. 


X 

The  units  of  watts  per  square  centimeter  and  joules  per  square 
centimeter  are  standard  for  laser  hazard  definitions  and  will 
be  used  here  and  in  Section  4.4.7, 
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Class  III:  For  a scanned  laser,  the  radiant  exposure 
through  a 7 mm  diameter  aperature  shall  not  exceed  10  joules/ 
cm  . For  a fixed  beam  laser,  the  radiant  power  through  an 
80  mm  diameter  aperture  shall  not  exceed  0.5  watt. 

Class  IV:  Does  not  meet  Class  III  requirements. 

In  its  intended  operational  modes,  the  laser  beam  will  always 
be  scanning.  Using  the  worst  case  of  a focal  volume  at  a 50 
foot  (15  m)  range.  Figure  4-11  shows  the  Class  I limits  and  the 
energy  levels  attained  by  a scanning  CW  laser  for  different 
emitted  powers.  Since  even  at  20  watt  laser  actually  emits 
less  than  10  watts  of  power,  due  to  transmission  losses,  it 
appears  that  CW  lasers  up  to  70  watts  can  be  utilized  without 
violating  Class  I standards.  (Class  III  standards  are  about 
3 orders  of  magnitude  more  lenient.) 

The  laser  system  should  have  provisions  to  shut  off  the 
transmitter  if  the  scanner  stops  during  operation.  With  these 
provisions.  Class  I limits  will  not  be  exceeded. 

4.4  PULSFD  LASERS  IN  THE  HORIZON  SCAN  MODE 
4,4.1  General  Description 

A pulse  transmitted  from  a laser  will  travel  along  a narrow 
beam  at  the  speed  of  light,  illuminating  at  every  instant  a 
tubular  volume  whose  length  is  proportional  to  the  pulse  width. 

The  received  signal  at  an  instant  some  time  after  the  transmission 
represents  the  sum  of  the  scattered  energies  from  that  volume 
illuminated  at  half  that  time  after  transmission.  By  examining 
the  return  signal  during  periods  of  duration  equal  to  the  pulse 
width,  a series  of  velocities  can  be  generated  for  each  transmitted 
pulse,  each  associated  with  a different  range  interval,  or  bin. 

This  is  discussed  in  detail  in  Section  4.4.5.  For  example,  a 
two  microsecond  pulse  is  associated  with  a 1000  foot  (300  m) 
range  resolution  cell.  By  examining  the  return  signal  during 
successive  two-microsecond  periods,  the  Doppler  velocities  can 
be  extracted  from  a series  of  1000  foot  (300  m)  range  bins.  Figure 
4-12  shows  a number  of  range  bins  and  the  relative  weighting 
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attached  to  signal  returns  from  particles  at  different  ranges; 
e.g.,  the  particles  at  19.5  kilofeet  range  contribute  ciiua  1 I y 
to  the  signals  in  bins  19  and  20,  and  not  at  all  to  bins  19  and 
21.  Thus  there  is  some  overlap,  but  75  percent  of  the  energy  in 
the  return  signal  is  from  within  +500  feet  (150  m)  of  the  range 
bin  center. 

For  a minimum  range  of  0.5  nm  (1  km)  and  a maximum  range  of 
7 nm  (13  km),  there  are  39  range  bins,  and  a doppler  velocity 
can  be  estimated  for  each  range  bin  on  each  transmitted  pulse. 

The  transmitted  beam  has  the  property  of  maintaining  its 
diameter  for  a considerable  distance;  for  example,  a 12  inch 
(0.3m)  diameter  telescope  generates  a collimated  beam  about  12 
inches  in  diameter  which  maintains  its  diameter  out  to  about 
2 nm  (4  km).  Beyond  that  it  gradually  spreads  to  about  3 feet 
(0.9m)  at  7 nm  (12  km)  range.  Since  the  received  signal  intensity 
is  inversely  proportional  to  the  beam  area,  the  resulting  loss 
of  signal  due  to  beamspread  is  about  9 dB  at  7 nm  (12  km)  range. 

Shears  will  be  evidenced  by  variations  in  the  Doppler  velocit- 
ies in  successive  ranges.  The  gradients  in  these  1 ine-of -s ight 
wind  velocities  are  essentially  the  headwind  shears  that  would  be 
experienced  by  an  approaching  or  departing  aircraft.  This  is 
due  to  the  geometry,  whereby  a sensor  located  near  the  center  of 
an  airport  has  radials  which  are  nearly  parallel  to  the  approach 
and  departure  lanes  (see  Figure  5-1). 

4.4.2  Range  Capability 

A pulsed  laser  system  will  not  achieve  the  desired  system 
range  of  7 nm  under  all  conditions:  under  conditions  of  heavy 
fog  or  moderate-to-heavy  rains,  the  range  will  drop  due  to  severe 
attenuation.  This  limitation  poses  the  largest  problem  to  the 
deployment  of  pulsed  laser  sensors.  Fortunately,  heavy  fogs  are 
usually  associated  with  warm  front  shears,  rather  than  thunder- 
storms; in  the  latter  case  some  inconsequential  light  fogs  may 
accompany  rain.  It  is  not  anticipated  that  pulsed  lasers  would 
be  used  to  detect  warm  front  shears. 
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The  features  of  thunderstorm  related  shears  that  are 
encouraging  are  that  (1)  gust  fronts  and  secondary  surges  of 
the  cold  air  outflow  precede  the  rainfall  of  a parent  thunder- 
storm, and  thus  can  be  readily  detected  in  advance  (see  Goff, 

1976);  and  (2)  downburst  cells  usually  have  heavy  rain  within 
the  cell  but  are  surrounded  by  areas  of  light  rainfall  or  no 
rainfall  (Fujita,  1975,  and  NTSB,  1974b,  1976a,  and  1976c). 

In  these  conditions  pulsed  lasers  would  be  able  to  achieve 
adequate  system  ranges.  The  evidence  is  not  as  extensive  as 
one  would  like,  and  needs  further  experimentation. 

The  effect  of  attenuation  on  range  can  be  seen  from  Figures 
4-13  and  4-14,  which  are  calculated  using  system  parameters 
which  arc  believed  to  be  readily  achievable  with  today's 
technology  (see  Appendix  A).  It  can  be  seen  that  a system  range 
7 nm  (13  km)  can  be  achieved  for  attenuations  less  than  2 dB/nm 
(1.1  dB/km)  (assuming  a SNR  threshold  of  3 dB) . Comparing  this 
to  Figure  4-4,  it  can  be  seen  that  this  corresponds  to  about 
2 mm/hr.  Thus  if  the  average  rainfall  along  the  beam  path  does 
not  exceed  2 mm/hr,  the  system  range  can  be  achieved  with  10 
kilowatts  of  transmitter  power. 

Higher  laser  performance  would  be  helpful  in  improving  the 
range  but  not  as  much  as  one  might  expect.  For  example,  a 1000 
percent  increase  in  transmitter  power  would  allow  system  ranges 
of  7 nm  in  average  rainfalls  of  up  to  3mm/hr  (0.1  in/hr). 

Another  1000  percent  increase  in  transmitter  power  would  increase 
this  to  about  4 mm/hr. 

An  interesting  point  to  note  is  that  a line-of-sight  sensor 
l)ointing  up  the  glide  path  will  not  detect  downburst  cells  beyond 
about  1-2  miles  from  touchdown  even  if  system  range  were  unlimited. 
This  can  he  seen  from  Figure  4-15:  above  the  outwash  of  a down- 
burst  cell  the  wind  currents  are  essentially  downdrafts,  which 
couple  weakly  into  the  line-of-sight  velocity  component.  This, 
of  course,  can  he  alleviated  by  pointing  the  laser  below  the 
glide  path.  It  is  suggested  that  the  beam  be  pointed  at  such 
an  elevation  angle  that  it  crosses  the  glidepath  at  about  250-300 
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feet  (75-90  m)  in  altitude,  as  shown  in  Figure  4-15.  It  can  be 
seen  that  the  outwash  of  the  cell  will  be  detected  by  the  laser 
sensor  at  much  longer  ranges  when  it  is  pointed  at  lower  angles. 


4.4.5  Scanning  Capability 

Before  scan  speed  and  data  update  rate  are  discussed,  the 
need  for  scanning  must  be  established.  First  of  all,  without 
scanning,  a sensor  would  have  no  advance  warning  capability.  It 
would  answer  the  question,  "Is  there  a wind  shear  event  in  the 
corridor?",  but  not,  "Is  there  a wind  shear  event  approaching  the 
corridor?"  Secondly,  unless  a separate  sensor  were  placed  at 
each  threshold,  a scanning  mechanism  would  still  be  required 
for  pointing  successively  at  different  approach  corridors.  To 
require  a separate  sensor  for  each  approach/departure  zone  would 
be  prohibitively  expensive  for  most  airports. 

For  wide-scale  synoptic  phenomena  such  as  inversions  and 
warm  fronts,  it  is  less  critical  to  have  advance  information. 

In  the  case  of  rapidly  advancing  cold  fronts  and  thunderstorm 
outflows,  however,  a hazardous  situation  could  arise  in  the 
two  minutes  it  takes  a typical  aircraft  to  travel  from  the  outer 
marker  to  decision  height.  In  these  two  minutes  a thunderstorm 
system  moving  at  30  knots  would  travel  a nautical  mile.  Since 
the  scan  capability  is  required,  anyhow,  it  should  be  used  to 
detect  the  incoming  hazard. 

The  discussion  on  update  rate  in  Section  3.2  concludes 
that  an  update  rate  of  two  per  minute  is  adequate  to  detect 
fast-moving  fronts,  gust  fronts,  and  downburst  cells.  Since 
the  storm  system  (or  front)  can  approach  from  any  direction,  it 
is  desirable  to  have  360  degree  coverage,  especially  until  a 
storm  system  and  its  motion  have  been  established.  The  scan 
rate  for  a two-per-minute  update  rate  is  360  degrees  in  30 
seconds,  or  12  degrees  per  second.  It  turns  out  that  speed 
of  light  limitations  do  not  allow  faster  scan  rate  if  the 
azimuti)  scanner  is  smoothly  rotated;  this  can  be  demonstrated  by 
the  following  argument:  the  half-power  beamwidth  of  a 
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laser  using  1 foot  optics  is  given  by  1.3  X/D  radians,  or 
0.046  milliradian.  At  a range  of  7 nm  (13  km),  it  takes  84 
microseconds  for  a round  trip;  thus  at  a scan  rate  of  12  degrees 
per  second,  the  scanner  has  moved  12x84x10*^-10’^  degrees,  or 
0.0175  milliradian,  or  3/8  of  a beamwidth.  This  represents 
a signal  loss  of  about  1 dB,  which  is  acceptable.  A telescope 
twice  as  large  would  lose  about  6 dB  of  signal  at  7 nm  (13  km), 
or  would  be  limited  to  a range  of  3.5  nm  (6.5km).  Thus  larger 
apertures  would  not  improve  system  performance  in  a scanning  mode. 

Fortunately,  the  fact  that  convenient  optics  sizes  are 
concomitant  with  the  desired  update  rate  means  that  a constant 
rate  of  scan  in  azimuth  can  be  employed  (otherwise,  a step  scan 
mode  would  have  to  be  used).  For  example,  a 12  degree  per  second 
scan  with  a laser  having  a pulse  repetition  rate  of  24  per  second 
yields  a mapping  with  0.5  degree  resolution,  which  corresponds  to 
370  feet  (113  m)  at  7 nm  (13  km)  range.  These  parameters  are 
quite  reasonable  for  the  application. 

4.4.4  Pulse  Width 

The  pulse  width  is  determined  by  a tradeoff  between  range 
resolution  and  velocity  resolution.  Figure  4-12  shows  the  rela- 
tive signal  weighting  attached  to  returns  from  particles  associ- 
ated with  a given  range  for  a two  microsecond  wide  square  pulse. 
While  returns  from  a 2000  foot  (600  m)  volume  figure  in  the  'total, 
75  percent  of  the  return  comes  from  the  central  1000  feet  (300  m) 
(see  Section  4.4.5  for  a full  discussion  of  this);  for  this  sensor, 
a two  microsecond  wide  pulse  is  said  to  be  associated  with  a 
1000  foot  (600  m)  range  resolution  cell. 

Velocity  resolution  is  directly  related  to  the  frequency 
resolution  of  the  return  signal  and  the  signal  processor. 

Generally  the  signal  processor  bandwidth  is  matched  to  the  pulse 
spectral  width.  For  a two  microsecond  pulse,  the  processor 
bandwidth  would  be  500  kHz,  and  for  a four  microsecond  pulse, 

250  kllz.  Roughly  speaking,  one  knot  of  speed  corresponds  to 
100  kHz  (actually,  97.16  kHz).  Thus  a 500  kHz  bandwidth 
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corresponds  to  about  S knots  of  resolution.  The  tradeoff  is 
shown  in  Figure  4-16. 

The  velocity  accuracy,  as  opposed  to  resolution,  depends  on 
the  signal - to-noise  ratio  as  well  as  the  signal  processor 

,V  = = 10^  knots 

2/SM  /SNR 

where  B is  the  bandwidth  in  megahertz.  For  example,  for  signal- 
to-noise  ratios  greater  than  14  dB,  the  theoretical  accuracy 
(per  pulse)  is  better  than  1 knot.  The  next  section  discusses 
the  implications  to  the  signal  processing  techniques  to  be 
employed.  The  important  conclusion  is  that  there  appear  to  be 
methods  available  to  achieve  wind  speed  accuracies  of  1-2  knots 
(0.5-1  m/sec)  with  two-microsecond  pulses. 

4.4.5  Averaging  and  Signal  Processing 

It  was  mentioned  in  the  previous  section  that  the  inherent 
accuracy  of  the  velocity  estimate  is  better  than  the  matched 
filter  resolution.  It  is  also  true  that  proper  averaging  of 
velocity  estimates  from  successive  pulses  will  improve  the 
final  velocity  estimate  for  the  range  bins  associated  with  a 
particular  azimuth.  The  purpose  of  averaging  and  signal  process- 
ing is  to  obtain  the  best  estimate  of  the  velocity  characteristic 
of  the  volume  of  interest. 

In  principle  a complete  velocity  spectrum  is  received  at 
every  range  for  each  transmitted  pulse:  at  every  instant  of 
time  after  the  pulse  begins  to  leave  the  transmitter  it  illuminates 
some  volume  along  the  beam  line  of  sight.  Assuming  a simple 
square  pulse,  the  volume  is  approximately  a cylindrical  tube, 
having  a beam  diameter  equal  to  the  telescope  diameter  initially, 
and  a length  equal  to  the  speed  of  light  times  the  pulse  width. 

I'or  a 2 microsecond  pulse  width,  the  tube  is  about  2000  feet 
(600  m)  long  at  any  instant.  The  situation  is  shown  in  Figure 
4-17a,  The  returns  at  tg,  for  example,  are  made  up  of  the  returns 
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FIGURE  4-17a.  TIME-RANGE  DIAGRAM 
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FIGURE  4-17C.  RELATIVE  WEIGHTING 
OF  SPATIAL  CONTRIBUTIONS  TO  PULSE 
RECEPTION  (RANGE  GATE  DURATION 
EQUAL  TO  PULSE  WIDTH) 
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from  Xj^  to  X2.  If  the  range  gate  is  left  on  for  a duration  longer  1 

than  the  pulse  width,  the  received  signal  ic  made  up  of  contribu-  1 

tions  from  X^  to  X^  weighted  as  shown  in  Figure  4-17b.  If  the 

range  gate  is  matched  to  the  pulse  (At  = t),  the  triangular 

weighting  function  of  Figure  4-17c  applies.  Here  75  percent  of 

the  return  signal  originates  from  the  central  portion,  ct/2  wide; 

thus  for  a 2 microsecond  pulse,  75  percent  of  the  energy  received 

during  the  range  sample  emanates  from  the  central  1000  feet  of  the 

volume. 

The  Doppler  shift  ranges  from  about  100  kHz  for  1 knot  winds 
to  6 Mllz  for  60  knot  winds.  It  is  anticipated  that  the  received 
signal  will  be  heterodyned  with  a local  oscillator  signal  at  a 
frequency  higher  (or  lower)  than  the  transmitted  frequency  in 
order  to  distinguish  positive  and  negative  velocities.  Assuming 

■ 

a frequency  offset  of  30  MHz,  there  will  between  24x2=48  and  I 

36x2=72  cycles  received  for  each  range  within  a 2 microsecond  j 

pulse.  This  is  a long  enough  record  to  enable  a spectrum  analy- 
sis for  each  pulse  and  each  range  bin. 

A filter  bank  having  filter  bandwidths  matched  to  the  pulse 
width  will  produce  at  any  given  time  a set  of  outputs  whose 
amplitudes  provide  a sampled  measure  of  the  velocity  spectrum  of 
the  corresponding  range  bin.  For  example,  a filter  bank  matched 
to  a two  microsecond  pulse  width  would  have  filters  500  kHz  wide; 

24  filters  would  cover  the  velocity  range  of  +60  knots  (+30  m/sec). 

If  40  range  bins  were  utilized,  the  filter  outputs  would  be  sampled 
40  times  at  2 microsecond  intervals.  The  same  considerations 
apply  here  as  for  the  CW  laser  (see  Section  4.3.6)  except  that  the 
velocity  resolution  here  would  be  almost  5 knots  (2.5  ra/sec) . 

Turbulence  and  wind  gradients  cause  a spreading  of  the  spectrum. 

For  example,  a significant  gradient  of  2.5  knots  per  thousand 
feet  (0.4  m/sec/km)  (see  Section  3.4)  would  have  a spectral  spread 
of  about  250  kHz,  wide  enough  so  that  significant  amplitudes  would 
be  found  in  two  or  three  adjacent  velocity  samples.  This  would 
allow  interpolation  to  estimate  the  velocity  centroid  to  better 
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than  the  nearest  5 knot  increment.  Other  techniques  could  also 
be  employed  to  this  data  granularity  problem. 


H 


n 

I 


One  approach  to  improving  the  accuracy  would  be  to  store  the 
waveform  from  pulse  to  pulse  and  analyze  the  extended  record. 

This  would  be  quite  beneficial  if  pulse-to-pulse  transmitter 
coherence  were  maintained,  and  if  the  returns  maintained  their 
phase  correlation  from  pulse  to  pulse.  However,  in  the  inter- 
pulse period  the  particle  configuration  changes  enough  to  destroy 
the  phase  correlation:  a mild  horizontal  shear  of  1 knot  per  1000 
feet  (0.16  m/sec/km)  results  in  a relative  motion  of  50  wavelengths 
per  millisecond;  thus  decorrelation  takes  place  in  microseconds, 
and  coherent  averaging  is  not  feasible.  Furthermore,  for  pulse 
repetition  rates  less  than  1000  pps,  the  sample  volume  has  com- 
pletely changed,  since  the  scanner  is  moving  in  azimuth. 


Even  with  a lack  of  phase  coherence,  averaging  is  beneficial. 
It  reduces  the  noise  variance,  thus  allowing  tighter  thresholding 
to  be  employed.  Incoherent  averaging  can  take  a variety  of  forms, 
one  of  which  is  to  store  the  amplitude  of  each  velocity  bin  (for 
each  range)  after  spectral  processing,  and  add  the  contents  from 
pulse  to  pulse  for  a period  of  time  determined  by  the  azimuth 
sampling  rate. 


Signal  processing  must  be  sophisticated  enough  to  reject 
spurious  measurements  and  measurements  whose  SNR  falls  below 
some  threshold,  e.g.,  3dB.  The  latter  situation  will  routinely 
occur  due  to  backscattering  variations,  described  in  Section 
4.2.3.  Other  improvements  will  be  possible  when  more  field  data 
become  available. 


Considering  the  benefits  of  averaging,  it  is  evident  that 
the  pulse  repetition  rate  should  be  as  high  as  laser  technology, 
budget  considerations,  processing  time,  and  safety  considerations 
will  allow;  also,  full  use  should  be  made  of  the  information 
available  on  each  pulse.  The  NASA  Clear  Air  Turbulence  (CAT) 
detector  employs  a laser  amplifier  operating  at  140  pulses  per 
second  (pps)  (Sonnenschein,  et  al,  1970).  A modern  off-the-shelf 


TEA  laser  typically  provides  2 pps,  while  20  pps  have  been  readily 
obtained  in  laboratory  setups;  although  there  are  technical  problems 
encountered  at  higher  repetition  rates  for  TEA  lasers,  there 
appears  to  be  no  fundamental  reason  why  200-1000  pps  cannot  be 
obtained.  Electron- injection  techniques  have  been  demonstrated 
which  achieve  megawatt  outputs  and  kilohertz  repetition  rates 
.Simultaneously . 

4.4.6  Data  Processing 

Data  processing  treats  the  samples  that  emerge  from  the 
signal  processor  and  is  concerned  with  the  reasonableness  of 
the  data  sample,  its  internal  structure,  and  the  sample-to-sample 
consistency.  Problems  that  data  processing  should  be  able  to 
handle  include  the  following: 

a.  Spectral  broadening  due  to  turbulence  and  shear. 

b.  Double-peaking  caused  by  heavy  objects  such  as  hail- 
stones, large  raindrops,  insects,  and  birds  that  are 
in  motion  with  respect  to  the  local  air  mass. 

c.  Spurious  errors  that  cause  data  outliers. 

Data  processing  also  includes  signature  recognition  algo- 
rithms which  involve  a priori  knowledge  of  wind  shear  descriptors, 
so  that  shears  can  he  distinguished  from  short  term,  large  scale 
gusts,  and  so  that'  downdraft  cells  can  be  distinguished  from 
frontal  shears.  It  is  difficult  to  state  precisely  the  require- 
ments of  such  algorithms  without  field  data.  Some  reasonable  map 
signatures  are  shown  in  Figure  4-18  for  some  well-defined  situa- 
tions, but  these  are  qualitative  onlv.  It  nav  be,  for  example, 
that  downdraft  identification  may  require  different  thresh- 
olds than  other  shear  types.  Turbulence  may  prove  to  be  highly 
correlated  with  certain  shear  regimes,  so  that  t*'*;  spectral  width 
variations  may  be  useful.  Certainly  the  amplitut  . the  return 
signal  is  related  to  the  precipitation  level  through  the  attenua- 
tion and  backscattering  mechanisms.  For  example,  a heavy  precipi- 
tation would  cause  the  signal  to  decrease  with  range  at  a much 
more  rapid  rate  than  would  be  true  in  light  rainfall.  This  in- 
formation would  be  used  to  help  identify  downdraft  cells,  which 


(♦)  INCREASING  HEADWIND  TO  APPROACHING  AIRCRAll 


(-)  DECREASING  HEADWIND  TO  APPROACHING  AIRCRAIT 


INVERSION  OR  WARM  FRONT 


b . COLD  FRONT 


THUNDERSTORM  GUST  FRONT 


d.  DOWNBURST  CELLS 


FIGURE  4-18.  PULSED  LASER  MAP  SIGNATURES 
OF  WIND  SHEAR  EVENTS 


usually  have  heavy  precipitation  within. 

Signature  recognition  is  clearly  an  area  where  more  work 
needs  to  be  done.  Field  data  obtained  during  shear  events  will 
prove  invaluable  here. 

4.4.7  Laser  Hazard  Limitations 

The  safety  requirements  for  a pulsed  laser  must  be  calcu- 
lated in  a manner  different  from  those  for  the  CW  laser.  For  one 
thing,  the  pulsed  laser  is  only  on  for  a burst  of  a few  micro- 
seconds: This,  coupled  with  the  fact  that  laser  damage  to  human 
and  animal  life  is  a function  of  average  power  rather  than  peak- 
power,  allows  quite  high  peak  power  levels  to  be  reached  safely. 

For  another,  the  pulsed  laser  beam  is  collimated;  the  CW  beam 
on  the  other  hand,  is  focused  on  a small  region,  wherein  high  local 
intensities  can  be  reached. 

Reasonable  safety  precautions  will  eliminate  situations 
where  accidental  irradiation  could  occur: 

a.  The  space  through  which  the  laser  beams  travel  within 

the  housing  should  be  enclosed  to  prevent  accidental 
insertion  of  hands  or  other  members  into  the  active 
region;  removal  of  the  cover  by  unauthorized  personnel 
should  automatically  shut  down  the  laser. 

b.  The  scanner  should  not  be  located  anywhere  that 
people  are  likely  to  be  working;  it  should  be  placed 
above  such  an  area. 

c.  The  scanner  should  have  mechanical  limits  to  prevent  it 
from  pointing  down. 

d.  Software  should  not  allow  the  scanner  to  stop  in  one 
place  for  longer  than  some  prescribed  time,  e.g.,  0.3-1 
second,  in  an  operational  environment.  For  test 
purposes  this  would  not  be  necessary,  since  personnel 
operating  the  unit  during  tests  would  be  trained. 

e.  Access  to  the  scanner  should  cause  the  laser  to  shut 
down. 
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First,  the  limits  of  an  operational  (scanning)  laser  will 

be  addressed.  Within  about  2-3  nautical  miles  of  the  laser,  the 

telescope  has  collimated  the  radiati^'n  into  approximately  a tube, 

whose  diameter  is  approximately  that  of  the  telescope,  assumed 

here  to  be  1 foot.  In  this  region  intensity  is  highest,  and  the  ' | 

potential  hazard  is  highest;  it  is  conservatively  assumed  that  ■ 

all  of  the  radiated  energy  is  in  this  tubular  beam.  For  a ! 

scanning  beam,  an  aperture  of  7 millimeters  diameter  is  utilized  ! 

for  the  standard  (see  section  4.3.8  and  FDA  (1975)),  which  states  j 

that,  to  meet  Class  I limits,  the  radiant  energy  detectable  within  j 

an  aperture  of  7 millimeters  diameter  shall  not  exceed  4.4xl0"^ 

1/4 

t joules,  where  t is  the  exposure  time  in  seconds.  Note  that 

the  measurement  aperture  is  small  compared  with  the  beam  diameter, 

unlike  the  case  for  CW  operation,  where  the  reverse  is  true.  If 

the  pulse  repetition  frequency  is  small  enough,  the  7 mm  test 

aperture  will  be  illuminated  only  once  per  scan.  At  higher 

repetition  rates,  the  beam  will  have  moved  less  than  a beam  width 

between  pulses,  and  several  illuminations  will  occur.  Figure 

4-19  shows  the  limits  on  peak  power  and  pulse  repetition  rate  j 

imposed  by  Class  I limits.  The  flat  portion  results  from  the  ; 

limit  of  energy  per  pulse,  and  the  sloping  region  from  energy  per  ' | 

scan.  The  long-term  limit  is  less  stringent  and  is  not  shown.  j 

From  Figure  4-19,  it  can  be  seen  that  the  present  TEA  laser  | 

operation  is  at  the  Class  I limits  within  10  meters  of  the  j 

scanner  (actually,  transmission  losses  between  laser  and  telescope 
output,  typically  3-5  dB,  bring  the  operation  well  within  the  i 

limits).  The  maximum  anticipated  operation  will  likewise  he 
within  Class  I limits  when  transmission  lasers  are  included  for 
ranges  greater  than  100  m-.  It  is  certainly  desirable  to  meet 
Class  I requirements,  since  the  precautions  that  are  required 
are  minimal.  However,  if  it  proves  necessary,  these  can  be 
exceeded;  operation  will  still  remain  within  the  Class  III  limits. 

Extra  precautions  and  justification  would  then  be  required. 
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PEAK  POWER,  WATTS 


I 


The  previous  discussion  applies  to  an  operational  sensor. 

For  experimental  purposes,  Class  III  (or  even  Class  IV  limits) 

can  be  applied.  For  fixed  beam  operation,  e.g.,  for  calibration, 

2 

Class  III  limits  (10  joules/cm  ) would  be  exceeded  only  for  powers 

greater  than  100  KW  and/or  pulse  rates  higher  than  500  pps.  If 

this  proved  to  be  a problem,  fixed-beam  operation  could  be  limited 

to  safely  low  repetition  rates.  ' 

! 


5.  SYSTEM  CONCEPT  AND  RECOMMENDATIONS 


5.1  PROPOSED  SYSTEM  CONCEPT  AND  ALTERNATIVES 

5.1.1  Preferred  System  Concept 

The  system  concept  proposed  he'fe  is  that  of  a single  centrally 
located,  scanning  laser  sensor,  which  will  provide  approach  and 
departure  corridor  protection,  and  advance  warning  of  shear  events. 
The  laser  scanner  is  lightweight,  typically  75-100  lb.,  while 
the  laser  itself  is  compact,  with  processing  electronics  which 
can  be  located  atop  a building,  or  on  a small  shelter  on  the 
airport  surface  (see  Figure  5-1). 

The  laser  is  a hybrid  pulsed/CW  unit:  the  pulsed  amplifier 
can  be  bypassed  such  that  the  primary  laser  oscillator  can 
radiate  directly,  in  a continuous  wave  (CW)  Velocity-Azimuth- 
Display  (VAD)  mode,  which  employs  a conical  scan  at  elevation 
angles  of  60-75  degrees  above  the  horizon.  The  system  operates 
periodically  (e.g.,  once  every  15  minutes)  in  the  VAD  mode  to 
measure  wind  speed  and  direction  at  several  altitudes;  this 
information,  which  could  be  incorporated  in  the  ATIS  meteorologi- 
cal broadcasts,  serves  a dual  purpose;  (1)  to  detect  any  synoptic 
vertical  shears,  and  (2)  to  provide  a two-altitude  (e.g.  surface 
and  1000  feet  altitude)  wind  report  to  assist  the  pilot  in  choos- 
ing his  approach  speed. 

In  the  pulsed  mode,  the  unit  scans  360  degrees  of  azimuth 
every  30  seconds. 

The  elevation  angle  is  programmed  to  intersect  the  glide 
path  of  each  runway  at  about  300  feet,  (see  Figure  5-2).  Note 
that  this  improves  downburst  cell  detection  at  ranges  farther 
out  than  1 nm  from  touchdown  (see  Section  4.4.3). 

Staggering  is  employed  in  range  to  reduce  the  effects  of 
quantization.  On  successive  pulses,  the  range  starting  point 
is  changed  by  half  a resolution  cell;  e.g.,  the  old  sweeps 
might  employ  range  bins  centered  at  1000  feet,  2000  feet,  3000 
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feet,  etc.  This  technique  helps  better  locate  and  assess  the 
shear  event. 

The  system  relies  on  line-sight  velocity  components  to 
derive  the  necessary  information.  Since  approaches  and 
departures  are  nearly  radial  outside  of  a mile  of  the  airport, 
a centrally  located  1 ine-of-sight  sensor  will  detect  primarily 
headwind  changes,  which  are  more  important  than  crosswind  changes 
(Figure  5-1) . 

The  data  processor  identifies  shear  zones  and  estimates 
their  intensity  and  spatial  extent.  Vertical  shears  are  identi- 
fied by  height  AGL,  and  are  calculated  for  each  corridor  orienta- 
tion. Horizontal  shears  are  located,  tracked,  and  estimates  are 
made  of  the  time  of  arrival  of  the  shear  event  at  each  corridor 
to  be  protected.  This  information  is  made  known  to  the  control 
tower.  The  manner  in  which  the  information  is  to  be  displayed 
has  not  yet  been  determined. 

5.1.2  Alternative  1 Laser  System;  Cross-Beam  Velocities 

Studies  conducted  to  date  (Lee,  1976)  indicate  that  spatial 
filtering  techniques  might  enable  a measurement  of  the  wind 
moving  across  the  beam.  This  would  enable  downdrafts  to  be 
measured  directly,  as  opposed  to  being  inferred  by  outwash 
signatures  in  velocity  and/or  amplitude  gradients  associated 
with  precipitation.  In  the  lateral  dimension,  the  downdraft 
cell  signature  would  be  enhanced  as  well.  However,  the  signal- 
to-noise  ratio  and  system  configuration  required  appear  to 
result  in  severe  range  limitations;  also,  the  theoretical  basis 
is  not  as  securely  established  for  this  approach  as  it  is  for 
the  conventional  1 ine-of-sight  approach.  It  is,  however,  an 
area  that  warrants  further  study. 

5.1.3  Alternative  2 Laser  System;  Dedicated  Corridor  System 

An  expensive,  but  possible,  alternative  is  to  dedicate  a 
laser  system  to  one  corridor,  or  one  approach-departure  corridor 
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pair.  The  advantage  of  this  is  that  even  under  heavy  rainfall  or 
fog,  the  system  range  would  still  be  1-2  nautical  miles,  and 
would  provide  excellent  corridor  protection.  By  eliminating  the 
scan  requirement,  some  cost  could  be  eliminated,  but  this  is 
minimal;  it  would  be  better  to  reduce  the  coverage  and  increase 
the  update  rate,  thereby  improving  system  performance. 

This  alternative  appears  to  have  no  drawback  or  special 
considerations  beyond  the  increased  cost.  It  would  be  located 
near  threshold  and  operate  at  an  elevation  angle  equal  to  or 
slightly  below  the  glideslope.  It  would  "see"  the  same  headwinds 
as  the  aircraft. 

This  alternative  would  be  impractical  for  an  application  to 
wind  shear  only;  it  has,  however,  the  capability  of  measuring 
vortices,  by  merely  adding  software.  Thus  it  could  be  an  element 
in  a Wake  Vortex  Avoidance  System.  'Additionally,  a remote 
visibility  sensor  could  share  much  of  the  same  equipment,  since 
it  would  be  colocated.  A multifunction  system  like  this  could 
be  economical  for  some  airports. 

5.2  OPERATIONAL  AND  DISPLAY  CONSIDERATIONS 

It  would  be  inappropriate  in  this  analysis  to  design  the 
operational  interfaces,  since  these  should  be  determined  by  the 
ultimate  users  of  the  system,  the  pilots  and  controllers.  Fuxther- 
more,  the  laser  system  is  flexible  enough  that  the  laser  concept 
need  not  be  tied  in  with  a particular  operational  procedure. 

However,  there  are  some  well -accepted  guidelines  which  are  useful 
to  mention  here,  that  are  important  in  the  development  of  the 
ultimate  sensor.  These  are  enumerated  below; 

a.  System  startup  should  be  a simple  switch.  No  personnel 
should  be  required  to  tune  or  adjust  equipment  at  start- 
up, nor  should  any  gas  fill  or  refrigerant  supply  be 
necessary.  Start-up  should  be  possible  by  a remote  switch. 
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Once  per  week  maintenance  for  about  half  an  hour 
should  be  sufficient.  If  laser  heads  are  to  be 
exchanged,  realignment  should  be  minimal,  not 
requiring  special  training.  A software  diagnostic 
routine  should  be  available,  and  be  simply  initiated. 

Displays  to  the  controller  should  not  require  constant 
survei 1 1 ancc ■ The  display  should  be  silent  and  blank, 
except  when  shear  events  have  been  detected.  When 
activated,  the  controller's  attention  should  be 
flagged,  probably  by  an  aural  cue. 

5.3  REQUIREMENTS  SUMMARY 

5.3.1  CW  System  Performance  Requirements 

Height  Range  - From  50  to  1000  feet  above  mean  ground  level.  j 

Height  Resolution  - 50  feet  or  less,  for  heights  up  to  j 

500  feet;  100  feet  or  less,  for  heights  of  500  to  1000  feet.  ! 

j 

Wind-Speed  Range  - +60  knots,  headwind  and  crosswind;  | 

+(1-20)  knots,  vertical  draft.  j 

Measurement/Averaging  Time  - A complete  set  of  range/angle  ‘ j 

measurements  shall  take  place  in  40  seconds  or  less;  the  averaging 
time  shall  not  exceed  20  seconds;  if  exponential  averaging  is  I 

used,  the  time  constant  shall  not  exceed  8 seconds.  j 

1 

Wind  Speed  Accuracy  - +1.5  knots,  all  components;  noise  +1  j 

knot,  all  components.  It  is  assumed  that  these  measurements  are 
taken  at  50-foot  height  intervals  below  500  feet,  and  include  the 
effect  of  any  averaging  over  small  height  intervals  below  500 
feet,  as  well  as  over  small' r height  resolution  cells. 

Wind  Shear  Accuracy  - +1  kphf  for  heights  up  to  500  feet; 

♦2  kphf  for  heights  from  500  to  1000  feet.  It  is  assumed  that 
t;hese  measurements  are  taken  at  50-foot  intervals  below  500  feet. 
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Scanning  Speed  - The  scanner  shall  be  capable  of  conical  scan 

at  speeds  from  0.025  to  0.5  Hz,  at  an  elevation  angle  of  60 

degrees  above  the  horizon. 

Environmental  Conditions  - The  equipment  shall  be  entirely 
enclosed.  The  system  shall  meet  the  specifications  of  5.3.1  under 
the  following  conditions: 

a.  30  knot  surface  winds 

b.  Rainfall  rate:  25  mm/hr. 

c.  Snowfall  rate:  80  mm/hr. 

d.  Fog:  RVR  of  1200  feet 

e.  Temperature:  -10  to  110  degrees  F 

f.  Humidity:  10  to  100  percent  RH. 

RFI  Environment  - The  equipment  shall  operate  without  degra- 
dation in  a severe  airport  noise  environment. 

Personnel  Hazard  - The  radiation  intensity  shall  not  exceed 
the  safety  limits  of  the  standard  for  lasers  (4.3.8)  for  airline 
passengers,  crew,  airport  maintenance  personnel,  operating 
personnel,  or  other  personnel  in  the  area.  If  radiation  safety 
limits  are  exceeded  at  the  scanner,  a safety  switch  shall  be 
installed  on  any  ladder  to  prevent  accidental  exposure  to  person- 
nel attempting  to  approach  the  scanner.  Range  scan  control  limits 
shall  be  employed  to  prevent  focusing  of  dangerous  radiation  at 
close  range;  similarly,  elevation  control  limits  shall  prevent 
scanning  below  1 degree  above  the  horizon.  A convenient  means 
of  absorbing  the  laser  radiation  during  extended  standby  periods 
shall  be  provided. 

5.3.2  Pulsed  System  Performance  Requirements 

Slant  Range  - 0.3  to  7 nautical  miles. 

Slant  Range  Rcsolution/Accuracy  - 1000  feet,  at  all  ranges. 
Accuracy  better  than  1 percent  of  range. 

Wind  Speed  - +60  knots,  line  of  sight,  with  wind  sense 
determined  by  the  system. 
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Wind  Speed  Accuracy  - +1.5  knots  between  adjacent  range 
cells;  +3  knots  absolute,  at  all  ranges  up  to  5 nautical  miles. 

Measurement/Averaging  Time  - Variable,  to  be  determined. 

Pulse  Repetition  Frequency  - 20  pps  required;  design  goal 
200  pps. 

Pulse  Width  - Nominally  2 microseconds. 

Peak  Power  - 50  KW  minimum;  200  KW  design  goal. 

Focus  - Nominally  set  for  infinite  range  (collimated  beam) 

Plevation  Angle  - Variable,  from  at  least  1 to  6 degrees 
during  operation. 

Azimuth  Scan  Rate  - 0.008  Hz  to  0.05  Hz,  average. 

Environmental  Conditions  - The.  equipment  shall  be  enclosed 
In  addition  the  system  shall  meet  the  specifications  above  unde 
the  following  conditions. 

a.  30  knot  surface  winds. 

b.  Temperature:  10  to  110  degrees  F 

c.  Humidity:  10  to  100  percent 

RFI  Environment  - The  equipment  shall  operate  without 
degradation  in  a severe  airport  noise  environment. 

Personnel  Hazard  - As  per  paragraph  5.3.2,  except  that  fhe 
standards  described  in  4.4.7  shall  be  employed. 


APPENDIX  A.  RANGE  CALCULATIONS 

The  signal-to-noise  ratio  (SNR)  for  a pulsed  system  is 
found  by  (Sonnenschein  and  Horrigan,  1971): 


where  : is  the  quantum  efficiency  of  the  detector 

: is  the  peak  power  of  the  transmitter,  in  watts 

^ 2 3 

B(ti):  is  the  backscattering  coefficient,  in  m /sr/m 

D:  is  the  telescope  diameter,  in  meters 

AR:  is  the  range  resolution,  in  meters 

B:  is  the  filter  bandwidth,  in  Hertz 

hv:  is  the  photon  energy,  in  joules 

is  the  transmission  efficiency  of  the  optical 

system 

R:  is  the  range,  in  meters 

X:  is  the  wavelength,  in  meters 
a:  is  the  attenuation  coefficient,  in  nepers 

per  meter. 

As  an  example  which  uses  values  representative  of  laser  technology 
today,  the  following  numbers  are  used: 

rij  ■ 50  percent 
Ti^  ■ 5 percent 
P^  - 10  KW 
D • 0. 3 m 
AR  - 300  m 
B - SOO  KHz 

hv  ■ 1.88  X 10*^®  joules 

X • 10.6  X lO'^  m 
2 3 

6(ti)  ^ a/80Tr  m /sr/m  (See  Appendix  B)  . 
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In  the  ncJiT  field,  R<<  and  the  SNR  becomes: 


SNR  - 2,5  X 10^  a 


In  the  far  field,  R>>  and  the  SNR  becomes: 
SNR  = 1.1  X 10^^  a e‘^“'VR^- 


The  signal-to-noise  ratio  is  usually  expressed  in  decibels: 

SNR  (dB)  - 20  logjp  (SNR) 

Note  that  SNR  is  a ratio  of  intensities,  and  therefore  unitless, 

The  attenuation  coefficient  a is  expressed  here  in  nepers 
per  meter;  more  convenient  units  are  decibels  per  kilometer 
and  decibels  per  nautical  mile.  They  are  related  to  a by  the 
following  relations: 

a(dB/km)  » 4343  a (nepers/m) 
a(dB/nm)  = 8043  a (nepers/m)  . 

These  formulas  and  values  are  used  to  determine  the  curves 
of  Figures  4-13  and  4-14. 
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APPENDIX  B.  THE  RELATION  BETWEEN  BACKSCATTER  AND  ATTENUATION 
COEFFICIENTS  FOR  WATER  DROPLETS 


The  calculation  of  backscatter  and  attenuation  (also  called 
extinction)  coefficients  are  both  derived  from  the  classic 
problem  of  scattering  of  a plane  wave  by  a dielectric  sphere. 

This  seemingly  simple  problem  is  in  fact  a complex  one,  and  many 
books  and  thousands  of  articles  have  appeared  on  the  subject. 

A typical  curve  of  total  scattered  power  as  a function  of  drop 
size  is  shown  in  Figure  B-1;  the  behavior  of  the  backscattered 
component  is  similar,  but  is  much  more  pronounced  at  the  resonant 
peaks  and  valleys.  If  backscattering  and  total  scattering  tracked 
precisely,  the  coefficients  would  be  related  by  a constant.  In 
the  regime  where  the  drop  diameter  is  large  in  terms  of  wavelengths, 
the  total  scattering  cross  section  of  each  drop  is  approximately 
twice  the  projected  area  of  the  sphere,  regardless  of  dielectric 
constant.  The  resonances  cause  variations  around  this  value. 

Since  drop  distributions  are  not  uniform,  we  may  assume  that  the 
resonant  effects  will  average  out.  A heuristic  argument  is  now 
presented  relating  average  backscattering  and  attenuation 
coefficients  for  diameters  large  compared  to  a wavelength. 

The  method  is  to  begin  with  the  ratio  of  backscatter  coeffi- 
cient to  extinction  coefficient  for  water  at  visible  wavelengths, 
and  successively  account  for  absorption  and  index  of  refraction 
changes  at  10.6  micron  wav-length. 

At  visible  wavelengths,  it  is  generally  agreed  that  the 
average  value  of  this  ratio  (5)  is  given  by  (Lifsitz,  1974); 

. A iill  . 0^ 

* a 4Tr 

2 3 

where  Bf/n)  is  the  backscatter  coefficient  in  m /sr/m  , and  a is 
the  extinction  coefficient  in  nepers  per  meter.  A plot  of  the 
ratio  from  another  set  of  calculations  is  given  in  Figure  B-2, 
taken  from  Twomey  5 Howell  (1965);  it  can  be  seen  that  4Tit  is 
well  approximated  by  0.6. 
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SCATTERING  CROSS  SECTION,  SQUARE  WAVELENGTHS 


FIGURE  B-2.  RATIO  OF  BACKSCATTERING  TO 
ATTENUATION  COEFFICIENTS  AT  VISIBLE 
WAVELENGTHS  FOR  A TYPICAL  WATER  DROPLET 
(FROM  TWOMEY  AND  HOWELL,  1965) 
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The  index  of  refraction  of  water  at  visible  wavelengths 
(0.03  micron  in  particular)  is: 

n(0.63)  = 1.33  + j lO'^, 
while  at  (,‘0,  vvavelengths  110.  b microns)  it  is: 

nllO.O)  = 1.18  + i 0.08. 

Ihe  imaginary  term  is  a measure  of  alisorption.  At  visil)le  wave 
lengths  it  is  negligible,  while  at  10.6  micron  wavelength  it  is 
quite  significant.  Absorption  is  the  dominant  cause  of  attenua- 
tion for  droplets  less  than  10  microns  in  diameter  (and  of  course 
contributes  nothing  to  backscatter) . 

•Since  the  total  scatter  coefficient  is  independent  of  the 
index  of  refraction,  the  question  reduces  to  that  of  the  relative 
backscatter  coefficient  with  different  indices  of  refraction. 

The  model  to  be  used  is  to  approximate  a spherical  drop 
hv  a cylinder  with  two  flat  faces: 


O 2a  r 

2a(]  n ) n^  , 1-jO. 


An  incident  plane  wave  from  the  left  of  unit  intensity  would 

{"-■-O' 

reflect  Vi+1/  units  ot  intensity  times  na  , the  projected  area, 
and  transmit  the  rest^^^^ j ^ 2^ i nto  the  drop.  At  visible  wavc- 

( n-lV 

lengths,  none  would  be  absorbed  in  the  drop,  and  Vn+l/  would  be 
reflected  by  face  2 back  through  the  drop;  transmitted  back 
through  face  1 would  emerge  an  intensity  equal  to: 


4n 


,(n+l) 
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As  the  diameter  varies,  the  backscattered  reflections  go  in  and 
out  of  phase,  resulting  in  large  variations  in  backscattered 
signal.  (Since  the  total  scattering  does  not  vary  appreciably, 
the  energy  is  scattered  in  other  directions  when  the  backscatter 
is  low.)  On  the  average,  the  total  intensity  is  just  the  sum  of 
the  two  reflected  intensities.  For  n = 1.33,  the  intensity  is 
proportional  to: 


I 


0393 


At  10.6  microns,  there  is  no  second  reflection,  since  all  energy 
entering  the  droplet  is  absorbed.  Thus 


I (10.6) 


■ 


0.0068 


for  n = 1.18.  The  ratio  between  these  is 

1(10.6)  ^ 

Additionally,  the  total  attenuation  is  due  to  scattering  and 
absorption;  half  the  energy  is  absorbed  at  10.6  microns,  while 
none  is  absorbed  ^t  visible  wavelengths.  Thus  the  ratio  of 
backscatter  to  attenuation  at  10.6  microns  is  approximately 

^(10.6)  =1.0^.  (0.17)  = or 

^ “ FO^  m^/sr/m^ 
if  a is  expressed  in  dB/km, 


B = 9.2  X 10‘^  a (dB/km)  m^/sr/m^ 


or  approximately. 


.-6 


B s 10  " a (dB/km) 


m^/sr/m^  . 


I 
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A number  of  points  from  results  by  Rensch  and  Long  (1970)  were 

checked,  and  the  ratio  c was  found  to  vary  between 
I 0 12  0 05 

'^nr  represents  some  middle  value.  Even  though 
the  assumption  breaks  down  in  fog,  representative  values  of  ^ 
for  fog  also  yield  values  of  about  For  drop  sizes  below 

2 microns  in  diameter,  absorption  exceeds  scatter,  and  the  ratio 
rapidly  decreases  for  decreasing  drop  size.  Thus,  in  a fog,  the 
backscattcr  is  primarily  due  to  the  drops  greater  than  2 microns 
in  diameter,  but  smaller  drops  can  cause  attenuation  by  absorp- 


APPENDIX  C.  BOUNDARY-LAYER  SHEARS  (HALLOCK,  1976) 
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The  surface  boundary  layer  is  that  region  of  the  atomsphere 
where  the  shear  stress  is  practically  constant.  The  region 
extends  from  the  surface  to  a few  tens  of  meters  in  altitude  on 
clear  nights  with  light  winds  and  to  about  two  thousand  meters 
in  altitude  on  sunny  summer  afternoons.  In  the  surface  boundary 
layer,  only  the  surface  conditions,  the  stability  condition,  and 
altitude  affect  the  wind  and  turbulence  structure.  The  logarith- 
mic wind  law  as  derived  from  the  Monin-Obukov  similarity  theory 
and  experiments  for  neutral,  stable,  and  unstable  conditions  is, 
in  simplified  form,  given  by: 


u 


(u*/k) 


where  u is  the  wind  magnitude,  u*  the  surface  friction  velocity, 
k the  von  Karmon's  constant  (k*0.4),  z the  measurement  height,  z^ 
the  roughness  length,  and  L the  Monin-Obukov  stability  length. 

The  function  is  an  empirically  derived  universal  function  of 

z/L : 


■ » 

♦(!)  ■ 5-2  f: 

z/L 

'"(t)  ? n-(i-i8^)'^/^}d^ij 

z/L 


(neut  rail 
(stable) 


(unstable ) 


for  every  stable  condition,  the  logarithmic  law  is  not  valid;  the 
layers  of  the  atmosphere  become  disconnected  and  large  scale 
frontal  motions  define  the  wind  profile  shape. 

The  net  effect  of  turbulence  in  the  atmosphere  is  the  pro- 
duction of  a thoroughly  mixed  state  by  the  transport  of  fluid 
properties  down  thei r j’espective  gradients.  In  the  friction 
layer  of  the  atmosphere,  the  turbulence  operates  to  transport 
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1 

I 

I 

1 


I 

L 


momentum  from  the  higher  layers  of  the  atmosphere  downward  toward 
the  ground  layer  to  replace  the  momentum  lost  through  frictional 
dissipation.  If  the  mixing  process  were  not  present,  tlie  wind 
shear  near  the  ground  would  become  very  large,  being  governed 
only  by  surface  drag  and  the  kinematic  viscosity.  Consequently, 
for  any  given  atmospheric  state,  the  maximum  value  of  wind  shear 
may  be  expected  to  occur  just  prior  to  the  onset  of  the  mixing 
process . 

Low-level  wind  shears  result  from  the  interaction  of 
synoptic  scale  flow  with  the  surface  of  the  earth.  The  boundary 
layer  can  be  segmented  into  a constant  wind  direction  layer 
(where  the  winds  behave  as  in  the  equations  above)  and  a wind 
turning  layer.  The  constant  wind  direction  layer  extends  from 
the  surface  up  to  about  150  meters  and  wind  shear  results  from 
the  boundary  condition  that  the  steady-state  wind  speed  must  be 
zero  at  the  ground.  The  turning  layer  occurs  above  the  constant 
wind  direction  layer  and  exhibits  a marked  turning  of  the  steady- 
state  wind  vector  with  an  increase  in  altitude. 

Consider  first  the  constant  wind  direction  layer.  When  the 
temperature  lapse  rate  is  neutrally  stable,  the  mechanical 
production  of  turbulent  kinetic  energy  is  neither  augmented  by 
the  buoyant  jiroduction  of  turbulent  kinetic  energy  (as  in  the 
unstable  case)  nor  suppressed  by  thermal  stratification  (as  in 
the  stable  case).  Altliough  strong  wind  shears  accompany  neutral 
conditions,  the  shear  is  confined  to  below  about  10  motors. 
Unstable  conditions  are  synonymous  with  convective  mixing.  Thus, 
both  neutral  and  unstable  conditions  should  not  produce  wind 
shears  that  are  hazardous  to  aviation.  .Stable  conditions,  on 
the  other  hand,  can  lead  to  dangerous  shear  conditions.  Negative 
buoyant  forces  suppress  the  turbulence  and  decouple  the  layers 
allowing  them  to  slip  relative  to  each  other.  Large  mean  flow 
wind  shears  can  then  result.  With  very  stable  conditions  the 
flow  is  laminar,  very  little  mixing  occurs,  and  the  overriding 
air  remains  separated  from  the  underlying  air.  Large  changes  in 
wind  direction  (in  excess  of  45  degrees)  and/or  magnitude  (in 
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excess  of  8 knots)  are  possible  in  altitude  changes  of  only  tens 
of  meters. 


Innumerable  wind  profiles  are  possible  in  the  turninn  layer. 
Stable  conditions  produce  the  most  danserous  shears  in  the  turn- 
ing layer  with  typical  anp.le  chanpes  of  SO  to  50  deprees.  The 
stable  turning  layer  is  the  least  understood  layer  due  to  the 
decoupling  of  the  layers  as  a result  of  the  decreased  turbulence. 
During  the  night,  an  inversion  can  build  up  and  relax  the  con- 
straint imposed  on  the  wind  by  day-time  mixing.  As  a result, 
the  wind  at  the  top  of  the  inversion  accelerated,  becomes  super- 
geostrophic,  and  oscillates  inertially--  a low  level  nocturnal 
jet.  The  low-level  jet  is  a thin,  well-defined  region  of  high- 
speed air  at  a typical  altitude  of  about  500  meters,  with  no 
indication  of  such  a wind  near  the  ground.  The  jet  begins  to 
build  up  in  the  late  afternoon,  reaches  its  maximum  in  the  middle 
of  the  night,  and  decays  in  the  early  morning.  The  nocturnal 
1 ’ jets  can  produce  s ipn i f icant 1 v large  values  of  low-level  wind 

I shear.  At  the  peak  of  the  jet,  the  winds  in  its  core,  between 

' 250  and  700  meters  in  altitude,  can  attain  lietween  40  and  70 

I knots,  decreasing  to  10  or  20  knots  between  1000  and  1400  meters, 

; and  to  0 at  the  ground.  It  is  a local  phenomenon  but  can  be  1500 

j kilometers  long  and  70  to  700  kilometers  wide.  In  the  United 

i States,  the  occurrence  of  the  nocturnal  jet  is  most  pronounced 

^ in  the  Great  Plains.  Tower  measurements  have  recorded  low-level 

jet  shears  in  excess  of  12  meters/second/30  meters. 
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APPENDIX  D.  ANALYSIS  OF  WIND  SHEAR  ERRORS  AND  THEIR  EFFECTS  ON 
SYSTEM  PERFORMANCE 


If  the  wind  shear  criterion  could  be  simply  characterized 
by  two  states  (i.e.,  Severe  and  Not  Severe),  the  problem  could 
be  classically  stated  in  terms  of  False  Alarms  and  Misses,  where 
a False  Alarm  indicates  a "Severe"  indication  when  the  actual 
shear  is  "Not  Severe"  and  a "Miss"  indicates  a "Not  Severe" 
wlicn  the  actual  shear  is  "Severe."  If  more  states  are  desired 
(i.e.,  I.iftht,  Moderate,  Strong,  Severe),  the  characterization  is 
more  complicated. 

I 

Staying  with  the  two-state  criterion,  it  is  convenient  to 
use  a mixed  continuous/discrete  description  whereby  wind  shear 
is  characterized  by  whole  kphf;  where,  for  example,  5 kphf 
means  a shear  between  4.5  and  5.5  kphf  . Let  Wj  be  the  thresh- 
old value,  and  W the  shear; 

|W(  > W.j.  implies  a Severe  shear,  and 

|W|  < Wrj.  implies  a Non-Severe  shear, 

where  W.j,  is  half  integer  value  (e.g.,  4.5  kphf,  5.5,  etc.) 

The  probability  of  a False  Alarm,  Ppy^i  is; 

PpA  = p(|WmI>Wt|iWaI<Wt)  n)-i) 

where  is  the  measured  shear,  and  is  the  actual  shear.  For 

example,  if  W.j,  = 4.5  the  event  |W^|<W.j.  means  that  can  be  -4, 

-3,  -2,  -1,  0,  1,  2,  3,  or  4 kphf;  however,  in  general,  the 
probabilities  of  these  subevents  are  not  uniform,  but  depend  on 
the  probable  density  of  the  phenomenon,  PjCWa)’ 

The  ICAO  paper  on  Low  Level  Wind  Shear  and  Turbulence  (ICAO, 
l!)7())  suggests  the  following  values  for  the  probability  of  exceed- 
ing the  values  shown; 
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3 kphf  - 50  percent 

5 kphf  - 1~  percent 

8 kphf  - 2 percent 

10  kphf  - 0.4  percent 


where  the  percentage  values  refer  to: 


J 


)dW^  - PjlWol 


(D-2) 


A little  curve-fitting  yields  the  following  aoproximation : 


p (WJ  = 0.09  exp{-0.105 

s A A A 


(D-3) 


and  P^(W  ) = 4 exp{-0.105!w  . 

So  2 o 


We  take  P..(m)  here  to  mean: 


p^(m)  4|p^(|m-|.|)  - Pg(|M+i|) 


(D-4) 


Thus  for  anv  event  A: 


l’(A  l‘Vl<W^)  = - 


[ ] 

Z W.»n)p  (n) 


(n-5) 


I 

n»-  [W.J,] 


where  (x)  is  the  largest  integer  less  than  or  equal  to  x. 
A is  the  event  that  |W|^|>W.j,.  Assuming  Gaussian  errors 
with  zero  mean  and  a standard  deviation  of  a kphf: 


where  E{e}  = 0,  and  E{e*}  = a*. 


(D-6) 


f.* 


I 


For  example,  suppose  = 3 kphf.  A False  Alarm  occurs 


VJ^=4.5 


|W_|=-4.5  0^ 


r! 


whf-rc  I’' ' IV.j! -Iv,  ll\^  = n)  is  given  by  (U-IO),  |i^(n)  is  given  In 
( I ) I j , and 


" = /. 

n n=-|l\j.l 


Similarly,  the  probabilities  of  a Miss  is  given  by 


I’m  = i’U.VmI<Wt 


= in 


IM.  |W^,|.W.,  W^=m1  Pj.(m) 


(I) 


I’s 


-lW^I-1 

where  T.*  = I 
m ;n  = - oo 


+ H 

m=lW,j,l+l 


, and 
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This  gives  the  following  results  for  W,p  = 4.5: 


0 

’’PA 

Pm 

1 

.087 

.154 

2 

.189 

. 257 

2.8 

. 261 

.307 

4 

. .560 

.3  50  . 

Changing  the  threshold  to  W, 


7.5  gives: 


GLOSSARY 


Aerosols  - Particles  suspended  in  the  atmosphere:  dust  parti- 
cles, moisture  droplets,  soot.  Sometimes  foji  and 
rain  droplets  are  also  called  aerosols. 

A(!l.  - Above  moan  Ground  Level;  used  to  denote  hei.ulit  above 

a runway.  It  is  used  in  contradistinction  to  alti- 
tuue,  which  is  referenced  to  mean  sea  level. 

Alls  - Automated  Terminal  Information  Service;  a recorded 

weather  advisory  for  pilots,  broadcast  from 
selected  airports. 

Hackscattcr  coefficient  - A coefficient  which  denotes  the  relative  ^ 

strength  of  the  laser  signal  reflected  from  a * 

volume  of  aerosols  back  toward  tlie  laser.  j 

boundary  layer  - The  lowest  kilometer  of  the  atmosphere;  so  named  ; 

because  the  effect  of  the  ground  is  to  retard  tlie 
winds  near  the  ground.  | 

Clutter  - The  unwanted  radar  returns  from  vegetation  and 

buildings;  clutter  tends  to  mask  the  returns  from 
desired  targets. 

Collimated  beani  - A beam  of  radiation  that  remains  in  a column 

for  some  distance  from  the  telescope  used  to  shape  , 

the  beam. 

(:02  laser  - A laser  operating  at  one  or  more  of  several  wave- 
lengths very  close  to  10.59  microns,  characteristic 
of  carbon  dioxide  molecular  vibrational  resonances. 

(TV  - Continuous  wave,  used  to  describe  a signal  which 

remains  on  for  long  periods  of  time;  as  opposed 
to  pulsed  signals. 

Downhurst  cells  - Tightly  confined  downdrafts  which  occur  as 

secondary  surges  behind  the  gust  front  of  a thunder- 
storm outflow;  they  frequently  contain  rain  (see 
Sec t ion  2.3). 

Clidepath  - The  vertical  guidance  provided  by  an  ILS;  generally 
it  is  a plane  oriented  at  about  three  degrees 
‘I  elevation,  and  intersecting  the  runway  at  the  GPIP. 

The  guidance  is  not  usually  reliable  below  100-200 
feet  AGL. 
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CIMI 


1 


Oust  front 


ilotcrodync 


Hit  - 


r.lide  Path  Intercept  Point,  the  |)oint  on  t Ito 
runway  centerline  that  the  glidcpath  would  inter- 
cept if  the  straight-line  guidance  were  projected 
to  the  ground.  It  is  usually  located  about  lOOO 
feet  from  threshold. 

The  leading  edge  of  the  cold  air  outflow  from  a 
thunderstorm;  sometimes  used  to  denote  the  zone 
behind  the  leading  edge,  where  gusty  winds  and 
vertical  drafts  occur. 

A technique  for  receiving  returns  from  a target, 
wherein  the  detector  mixes  signals  from  the  return 
signal  with  a local  oscillator  signal  to  achieve 
high  detection  efficiency. 

Used  here  to  denote  the  detection  of  a shear  event 
by  a laser  sensor  during  one  scan.  Hits  on 
successive  scans  i>rovide  a track  of  a shear  event, 
and  allow  an  estimate  of  its  velocity. 


Horizontal  sliear  - In  its  ideal  form,  a wind  condition  wherein 
wind  velocities  have  gradients  in  the  horizontal 
plane,  liut  are  constant  in  a vertical  line  above 
tlie  ground;  used  here  to  denote  a condition  where 
horizontal  gradients  are  larger  than  vertical 
grad i cuts . 

1 1.S  - Instrument  l.anding  .System,  a guidance  system 

associated  with  a runway,  which  provides  deviation 
from  runway  centerline  and  glidepath  to  instrumented 
aircraft.  An  outer  marker  is  generally  employed 
about  5 miles  out,  which  provides  timing  informa- 
tion; a middle  marker  is  also  employed,  loc.-ted 
about  half  a mile  out. 

Isotacli  - Contours  on  a wind  map  showing  lines  of  constant 

wind  speed. 


KIAS  - Knots  of  Indicated  Air  Speed,  the  speed  of  the 

aircraft  relative  to  the  air  mass. 


kjihf  - Knots  Per  Hundred  Feet,  a unit  used  here  to  describe 

the  intensity  of  a vertical  shear;  2 kphf  is  mild, 

5 kphf  is  significant,  and  12  kphf  is  treacherous. 


Localizer  - The  transmitter  used  in  an  ILS  to  provide  lateral 
guidance  in  the  form  of  deviations  from  runway 
centerline;  it  is  generally  located  about  1000 
feet  beyond  the  stop  end  of  the  runway. 
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Matched  filter  - A filter  in  a signal  processo  r/ rece  i ver  wliich 
is  shaped  to  optimize  the  SNR  to  a known  trans- 
mitter pulse  shape. 

Middle  marker  - A transmitter  in  an  ILS  which  provides  a signal, 
generally  audible,  to  the  pilot  of  an  approaching 
aircraft  that  he  is  about  half  a nautical  mile 
from  the  runway. 

Nocturnal  jet  - A thin,  well-defined  region  of  high-speed  air 
that  occurs  at  about  1000  feet  AGL  in  the  late 
afternoon  and  night;  it  would  not  be  detected  liy 
an  anemometer  near  the  surface. 

Outlier  - A measurement  wliich  differs  so  much  from  otlier 

measurements  proximate  in  space  or  time  that  it  is 
highly  unlikel)'  to  be  valid,  and  is  probably 
spur i ous . 

Rll  Kel.itive  Humidity 

•Scan  - A complete  set  of  measurements;  in  the  case  of 

a lasei'  scanner,  a scan  denotes  a complete  cycle 
of  range  and  angle  variations. 

SNU  - .Si  gnal  - to- No  isc  Ratio,  a me.isure  of  the  goodness 

of  the  signal  received;  low  SNR's  result  in  loss 
of  accuracy  and  increased  uncertainty. 

Sweep  - A set  of  measurements  associated  with  a single 

pulse,  in  a pulsed  laser  or  rad.ir;  the  set  of 
sweeps  comprising  the  azimuth  scan  form  a scan. 

System  range  - Used  to  denote  the  maximum  distance  from  the 
sensor  system  for  which  a useful  signal  is 
obt.ained;  as  opposed  to  "range,"  which  is  the 
distance  from  the  sensor  to  some  point  implied 
by  the  context. 

VAl)  - Velocity  Azimuth  Display,  ;■  scanning  measurement 

technique  used  by  CW  radars  and  lasers  to  obtain 
wind  profiles  (see  Section  4 . .'S ) . 

Vertical  shear  - a wind  condition  wherein  horizontal  winds 
exhibit  significant  gradients  in  the  vertical 
direction,  but  change  gradually  horizontally. 
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